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Abstract 

We perform a thorough study of thermal leptogenesis adding finite temperature 
effects, RGE corrections, scatterings involving gauge bosons and by properly 
avoiding overcounting on-shell processes. Assuming hierarchical right-handed 
neutrinos with arbitrary abundancy, successful leptogenesis can be achieved 
if left-handed neutrinos are lighter than 0.15 eV and right-handed neutrinos 
heavier than 2 x 10 7 GeV (SM case, 3a C.L.). MSSM results are similar. 
Furthermore, we study how reheating after inflation affects thermal leptoge- 
nesis. Assuming that the inflaton reheats SM particles but not directly right- 
handed neutrinos, we derive the lower bound on the reheating temperature to 
be Trh^2 x 10 9 GeV. This bound conflicts with the cosmological gravitino 
bound present in supersymmetric theories. We study some scenarios that avoid 
this conflict: 'soft leptogenesis', leptogenesis in presence of a large right-handed 
(s) neutrino abundancy or of a sneutrino condensate. 



1 Introduction 



If some new physics violates lepton number C at an energy scale Ac, neutrinos get small 
Majorana masses via the dimension-5 effective operator (LH) 2 /Ac- Experiments suggest 
Ac ~ 10 14 GeV. Indeed the solar and atmospheric data can be explained by neutrino 
oscillations induced by the following neutrino masses and mixings P 

|Am 2 t J = (2.0l8:f) x 10- 3 eV 2 , sin 2 2fl atm = 1.00 ± 0.04, 
Am s 2 un = (7.2 ± 0.7) x 10~ 5 eV 2 , tan 2 6 sun = 0.44 ± 0.05. 

Experiments will make further progress towards measuring effects accessible at low energy, 
completely described by 9 Majorana parameters: 3 neutrino masses, 3 mixing angles, 3 
CP-violating phases. 

One possible mechanism to generate the dimension-5 operator (LH) 2 / Ac is known as 
'see-saw' mechanism Adding three right-handed neutrinos N\^,3 with heavy Majorana 
masses mjv 3 > WJjv 3 > Mz and Yukawa couplings Y% 

& = ^sm + {^N 2 + Y-LiNjH + h.c.) , (2) 

one obtains light neutrino states with Majorana masses m„ = — (vY u ) T m^ 1 (vY u ). The 
see-saw is a simple and elegant mechanism, but hard to test experimentally. It predicts 
no relation between the 9 low-energy parameters, just reproducing them in terms of 18 
high-energy ones. The right-handed neutrinos which constitute the essence of the see-saw 
are too heavy or too weakly coupled to be experimentally observed. There are few possible 
indirect probes. For instance, in some supersymmetric models the Yukawa couplings 
might induce sizable rates for lepton flavour violating processes such as \i — > &y. On the 
cosmological side, thermal leptogenesis [3] provides an attractive scenario for the generation 
of the baryon asymmetry of the universe [3] . The three necessary conditions for the genera- 
tion of the baryon asymmetry [5] are satisfied in the Standard Model (SM) with additional, 
heavy singlet right-handed neutrinos: the baryon number is violated by sphaleron processes 
which convert the lepton asymmetry induced by the Majorana nature of the right-handed 
neutrino masses into baryon asymmetry; CP- violation is due to the Yukawa interactions of 
the right-handed neutrinos with the SM lepton doublets and out-of-equilibrium is induced 
by the right-handed neutrino decays. For some recent analyses see [HI El 01 ITOj. 

The goal of this paper is to perform a thorough analysis of thermal leptogenesis within 
the SM and the Minimal Supersymmetric Standard Model (MSSM). We improve the com- 
putation of baryon asymmetry generated through the mechanism of thermal leptogenesis 
by 

i) including finite temperature corrections to propagators, masses, decay and scattering 
processes, and to the CP-asymmetry; 

ii) renormalizing couplings at the relevant scale (~ 2ttT, where T is the relevant tem- 
perature, rather than ~ Mz); 

in) adding AL = 1 scatterings involving gauge bosons which turn out to be comparable 
or larger than the ones involving the top quark included in previous computations; 
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iv) performing a proper subtraction of washout scatterings mediated by intermediate on- 
shell particles (once this is correctly done, they are no longer resonantly enhanced); 

v) extending the analysis to situations where right-handed (s)neutrinos give a sizable 
contribution to the total energy density; 

vi) discussing how the predictions of thermal leptogenesis depend upon the cosmological 
assumptions. In particular, we study the effects of reheating after inflation and 
compute the lowest value of the reheating temperature T RH for successful leptogenesis. 

The paper is organized as follows. In section 121 we briefly summarize some general re- 
sults of field theory at finite temperature. In section El we discuss how we include thermal 
corrections and how they affect the different ingredients of leptogenesis. Details can be 
found in a series of appendices: Boltzmann equations in appendix scattering rates in 
appendix lB| CP-asymmetries in appendices IU1 (SM) andlDl (MSSM). In section |U we com- 
bine all ingredients to get the final baryon asymmetry predicted within SM leptogenesis, 
and study which thermal corrections turn out to be numerically important. On the basis 
of the lesson learned for the SM, in section E] we address the more involved case of super- 
symmetric leptogenesis. We apply our improved computation also to the 'soft leptogenesis' 
scenario ^] . In section El we discuss how the baryon asymmetry changes when the 
maximal temperature reached by the universe after inflation is not much higher than m^ l . 
The variation depends on one extra parameter, the reheating temperature Trh- Finally 
our results are summarized in section 

2 Finite-temperature propagators 

In order to consider finite temperature effects in the plasma, we work in the so-called real 
time formalism (RTF) ^3] of thermal field theory. The Green's functions computed in this 
formalism are directly time-ordered ones. The RTF requires the introduction of a ghost 
field dual to each physical field leading to the doubling of degrees of freedom. The thermal 
propagator has therefore a 2 x 2 structure: the (11) component refers to the physical field, 
the (22) component to the corresponding ghost field, with the off-diagonal components (12) 
and (21) mixing them. 

Scalars 

The complete propagator of a scalar particle (e.g. the Higgs) in momentum space is 

S®)-'™^ (3) 

coshfltf = v 7 ! + /bH , smh6 K = y/ f B (w), !b = H/r _ p (5) 
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where K is the particle four-momentum. In a general frame where the thermal bath has 
four-velocity (U^U^ = 1), we define the Lorentz-invariant quantities 



u = k = yJiKWtf - K»K^ (6) 

These coincide with the particle energy and momentum in the rest frame of the thermal 
bath, U>* = (1,0,0,0). 

In Eq. (JHJ), Ab{K) is the resummed propagator 

where tub is the bare mass and T,b{K) is the finite-temperature self-energy of the 
scalar boson field. Y>b{K) describes the continuous interactions with the heat bath al- 
tering the propagation of the boson, and it modifies the dispersion relation, substitut- 
ing particles with quasiparticles. At one- loop the self-energy takes the form Ti B {K) = 
Re T,b{K) + ilmUg^K), where KeHsiK) = m 2 B (T) is the effective plasma mass squared 
and 1m.H B {K) = — 2uoTb is proportional to the damping rate of the boson in the 
plasma. Since |ImS B (i^)| is suppressed compared to \ReT, B (K)\ [Hj, in the following 
we will work with resummed propagators for scalar bosons neglecting the absorptive part. 
Notice also that since we are considering one-loop thermal corrections to processes where 
all external fields are physical, we only need the (11) component of the bosonic propagator 



G n (K) = cosh 2 ^ — — + sinh 2 ^ 



— ? 



K 2 - m B {T) 2 + iej \K 2 - m B {T) 2 - ie 

where we have included the bare mass in mg(T). Using the property 

P ( - j - m5(x) , (9) 



x + ie \x 

where P denotes the principal value, the propagator can be rewritten as: 

G l \K) = — \ , + 2nf B (uj)5 [K 2 - m B (T) 2 } . (10) 

Fermions 

The (one-loop) resummed propagators for fermion fields can be written in RTF as 

S(K) = ( f:f K ] ) = M (T, K) ( 4 f A*f(K) ) M ^ V ■ (") 



where 



M(T,K)=[ C ° S t K ~ Sm r ) (12) 

SUMpK COS(pK 
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and 



cos0x = [0(lj) - 9(-u)] a/1 - f F (u) , sin(/) K = a/TfM 
In Eq. (JTTJ) A F (K) is the resummed propagator 

A F (K) 



e M/ T + 1' 



(13) 



T ^ - m - Ej,(tf) 



(14) 



where m is the fermion bare mass and H F (K) is the self-energy of the boson field at finite 
temperature. At one-loop the fermionic self-energy is given by ^Hj 



E F (K) = -a{K) -fK„ - b(K) -fU, 



(15) 



where is the four-velocity of the thermal bath as seen from a general frame. Neglecting 
the zero-temperature mass m , the coefficients a(K) and b(K) are given by 



a(K) 
b{K) 



m F (T) 
k 2 

m%(T) 
k 



uj , uj — k ujk T k 2 T 



UJ 
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(c^ 2 - fc 2 
2^ 



(16) 
(17) 



o n 2 k 



% [/b(p) + /Hp)] 



u; , 4p(p + fc) + k 2 - lu 2 , Ap 2 - (k + uj) 2 
m ^ o+P m 



2 4p(p-fc) + F-cj 2 ^ 4p 2 - (fc-cj) 2 _ 
7 I™ d P \f ( \ ft \U 4 P(P + fe) + A: 2 - c^ 2 

Here /b^ are the Bose-Einstein and Fermi-Dirac distributions, respectively, and m F (T) 
is the effective mass of the fermion in the plasma. Notice that this mass is not a "true" 
fermionic mass since it does not affect the chiral symmetry. The thermal mass is generated 
radiatively and comes from terms which are chiral-symmetric in the Lagrangian. The 
integrals / and J are gauge-dependent and have been computed here in the Feynman 
gauge. In the high-temperature limit (T ^> k), the terms in eqs. (|16|) and f)17|) proportional 
to the integrals / and J can be neglected, as they only give subleading contributions, leaving 
a result for the coefficients a and b which is gauge independent 1 . 

For a fermion charged under an SU(iV) gauge group with coupling g and having Yukawa 
coupling Y, the thermal fermionic mass squared is m F (T) = g 2 T 2 C{R) /8 + N f \Y\ 2 T 2 /16, 
where C(R) is the quadratic Casimir of the fermionic representation {e.g. C(R) = (N 2 — 
1)/(2N), when R is a fundamental of SU(N)) and Nf is the particle multiplicity flowing 
in the loop. 

Interactions of the fermions with the thermal bath modify the fermionic dispersion 
relation JH] leading to two different types of excitations with positive energy: 'particles' and 

■"■Our subsequent computations involve weakly coupled particles, so that m 2 F /T 2 turns out to be small 
and neglecting I and J is an excellent approximation. 
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Figure 1: The dispersion relation \w{k)\ (fig.^fl) and the residue (fig-W?) of particle (dotted 
line) and 'hole' (dashed line) excitations of a fermion with thermal mass m at temperature 
T for m < T . The solid line shows the approximation u> 2 = m 2 + k 2 . 



'holes' with the wrong correlation between chirality and helicity in the bare massless limit. 
The names 'particles' and 'holes' are suggested by an analogy with superconductors [TKj . 
This is because the propagator in eq. ()14j) has poles at uo = ±k — b/ (1 + a), an equation with 
two different positive-energy solutions E p and (E = see fig. ^i). At low- momentum, 

k <C m F , 

Ep = mF + ^ + 3^ + '"' ^ = mF ~^ + 3^ + '"' (20) 

O Olllp O Olllp 

while at larger momenta T ^> k ^> mp 

E p = k + ^-^\n^ T + --- , E h = k + 2kexp(-l-2k 2 /m 2 F ) + --- . (21) 

rC ZiK 771 jp 

The most remarkable property of hole dispersion relation is that its minimum occurs at 
the nonzero momentum ~ 0.4 mp. The residues of the particle and hole propagators are 
given by [To] 

E 2 — k 2 



^-b-(l + a){k-u) 
duo 



2m 2 F 



(22) 



which differs from the standard value Z = 1. This reduces to Z p = Zh = 1/2 at k mp and 
Z p = 1 and Zh = at k 3> mp (see fig. H]p): holes interact only at low momentum. Notice 
that the dispersion relation and the residues of the holes can be directly obtained from those 
of the particles, Eh(k) = E p (—k) and Zhik) = Z p (—k). The functions a(K) = a{k,uj) and 
b(K) = b(k,u) defined in eq. |T6j) - |T7j) . when evaluated for on-shell particles or holes are 
functions, for instance, only of k. Correspondingly, a useful property is a(k,E p (k)) = 
a(-k,E h (-k)) and \b(k,E p (k))\ = \b(-k, E h (-k))\. 

Despite these subtleties, since Z p + Z^ ps 1 the fermionic dispersion relation can be 
well approximated by u = [\k\ 2 + m^(T)] 1 / 2 (see fig. UJ and this is the form we have used 



5 



in our calculation, verifying that it is accurate enough when its validity seems doubtful. 
Notice that temperature-dependent masses enter into the kinematics and into the dispersion 
relations, but spinor functions must be taken identical to vacuum spinors, although with a 
modified dispersion relation |15j . 

As for the bosonic case, we are interested only in the (11) component of the resummed 
fermionic propagator 

S U (K) = cos 2 A - A F (K) - sin 2 x A* F (K), (23) 



A (K)= (1 + aYfK, + b^U, , A) 

n ' (1 + a) 2 K 2 + 2(1 + a)bK ■ U + b 2 + ie' 1 ' 

Assuming that the self-energy £ is real (which is always true in the high-T limit), so that 
a and b are real, with the help of eq. (jHJ), we can rewrite eq. (}23j) as 



[(1 + a)u + b] 2 - (1 + a) 2 k 2 + ie 
-2nf F (u)5 ([(1 + a)u + b] 2 - (1 + afk 2 ) 



(25) 



Before applying these results to the computation of the baryon asymmetry, let us briefly dis- 
cuss the issue of infrared singularities at finite temperature ^H]- It is well-known that per- 
turbation theory at finite temperature is afflicted by infrared problems which are worse than 
the ones appearing at zero-temperature field theory where the Kinoshita-Lee-Nauenberg 
(KLN) theorem QH] demonstrates that singularities appearing at intermediate stages 
of the calculation cancel out in the final physical result. In an interacting scalar theory the 
plasma mass may receive large two-loop infrared contributions and an all-loop resummation 
is needed to get a final result. In this case the plasma mass obtained at one-loop is used 
as an infrared regulator. In order to deal with infrared divergences one can perform the 
so-called hard thermal loop (HTL) resummation jTHj. In the HTL resummation one makes 
a distinction between hard momenta of order T and soft momenta of order gT and per- 
forms a resummation only for the soft lines. Indeed, the corrections to the bare propagator 
K 2 , being of order g 2 T 2 , start to be relevant only when K ~ gT, which is the soft scale. 
Therefore, if the internal momentum is hard, then ordinary bare propagators and vertices 
are sufficient, but if the momentum is soft then effective propagators and vertices must 
be used. In this way, one obtains an improved perturbation theory, in which the HTLs 
become generic analytic functions of the external momenta, g n T 2 f(ui, kj). Using the HTL 
technique, we have for instance explicitly checked that the top Yukawa coupling constant At 
entering the vertex HQ3U3 (which is involved in the AL = 1 scattering LH — > Q3U3) gets 
renormalized at finite temperature by the exchange of gauge bosons between the Higgs and 
the top lines and that the correction is tiny, 8X t /Xt — lO^ 1 g^- Furthermore one can extend 
the KLN theorem at finite temperature, and the infrared divergences in the corrections 
from virtual gauge-boson exchange cancel once we include absorption and emission of real 
particles that are too soft to escape the thermal bath 
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Figure 2: Thermal masses in the SM (left) and in the MSSM (right), in units of the 
temperature T. 

3 Including thermal corrections 

Performing a complete study of leptogenesis including all finite-temperature corrections is a 
very difficult task. In our computation we include the leading finite-temperature corrections 
which are quantitatively relevant for the computation of the final baryon asymmetry. These 
are given by: i) thermal corrections to gauge and Yukawa couplings; ii) thermal correc- 
tions to lepton, quark, and Higgs propagators; Hi) thermal corrections to the CP- violating 
asymmetry. 

In this section we discuss how we implement these three kinds of corrections and how 
large they are. Explicit formulas can be found in the appendices. In the next section we indi- 
cate how the finite-temperature computation of the baryon asymmetry can be significantly 
simplified by including only those effects which a posteriori turn out to be numerically most 
relevant. 



3.1 Corrections to couplings 

Renormalization of gauge and Yukawa couplings in the thermal plasma has been extensively 
studied (see e.g. ref. |21j). A very good approximation is to renormalize the couplings at 
the first Matsubara mode, 

E r = 2vrT, (26) 

using the zero-temperature renormalization group equations (RGE). This result can be 
understood by recalling that the average particle energy in the thermal plasma is larger 
than the temperature, and so must be the renormalization scale. For our purposes, the 
important couplings are the gauge and top Yukawa, which we evaluate using the appropriate 
RGE at the scale in eq. (J26|) . Therefore those couplings are always functions of temperature, 
even if not explicitly indicated. 
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1 10 2 10 4 10 6 10 8 10 10 10 12 10 14 10 16 
RGE scale \i in GeV 

Figure 3: Universal running of m v in the SM and in the MSSM. The bands give an indi- 
cation of the uncertainties, as explained in the text. 

Leptogenesis also depends on neutrino couplings and masses. The neutrino mass ma- 
trix can be renormalized from low energy (where it is measured) up to the high-energy 
scale relevant for leptogenesis using the well known RGE reported in appendix El They 
are often solved with a 'diagonalize and run' approach [2] which focuses on the neutrino 
masses and mixings probed by oscillation experiments. We instead purse a 'run and diag- 
onalize' strategy, as e.g. in ref. jS], which makes it easier to see how the combinations of 
neutrino masses relevant for leptogenesis renormalize. The solution can be trivially written 
as m(jl') = r I ■ m(p,) ■ I where r is an overall rescaling factor and 

/ ~ H + diarf v 2 x { " 3/2 in the SM (27) 

diag(y e , y^, y T ) ^ xj 1/cos2/5 in the M SSM 

(higher powers of ln(/x'//i) can be easily resummed). Here ye = m^/v with v = 174 GeV. 
Unless one considers large values of tan /3 the flavour dependent term / is very close to the 
unity matrix. In such a case, I can be relevant only if one considers special neutrino mass 
matrices, fine-tuned such that 2 or 3 eigenvalues are almost degenerate. We can neglect /, 
as long as we do not consider such cases. 

The relevant correction to neutrino masses is therefore given by the overall rescaling 
factor r. Its numerical value is plotted in fig. El The SM value has been computed assuming 
a 3 (M z ) = 0.118 ± 0.003, a pole top mass of m t = 175 ± 5 GeV and a Higgs mass m h = 
115 GeV. The band indicates the present uncertainty induced by the errors on m t and a^. 
Varying the Higgs mass has a negligible impact, unless is close to the triviality bound, 
m/i ~ 180 GeV, where the quartic Higgs coupling becomes non-perturbative at high scales 
inducing arbitrarily large values of r. This issue was also discussed in [23]. 

The MSSM central value has been computed assuming also moderately large tan (3 ~ 10, 
unification of gauge couplings at Mqut = 2 x 10 16 GeV, and A^Mqut) = 0.6. Ref. 
explains why A^Mgut) is still significantly uncertain, about between 0.5 -j- 0.7, giving rise 
to a correspondingly large uncertainty in RGE effects, illustrated by the shaded area in 

fig. El 
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3.2 Thermal corrections to decay and scattering processes 

To calculate the generated baryon asymmetry of the universe in the thermal leptogenesis 
scenario we have to solve a set of Boltzmann equations, discussed in appendix EJ which 
take into account processes that create or wash-out the asymmetry. We have recalculated 
the relevant reaction densities taking into account the propagation of particles inside the 
thermal plasma. 

As discussed in section the wave-functions that describe external fermion states in 
decay and scattering amplitudes are the same as in the zero-temperature case [IS]. Tem- 
perature corrections appear only in internal fermion lines and in the kinematics. In the 
SM, for all processes except those involving gauge bosons (which will be discussed later), 
the only fermion which mediates an interaction relevant for leptogenesis is N%. Since N% 
has no gauge interactions and its Yukawa coupling is small in all the relevant region of 
parameters, thermal corrections to Ni propagation can be neglected 2 . 

Therefore we only have to include thermal corrections to the dispersion relations of lep- 
ton doublets, third-generation quarks and Higgs bosons (and their supersymmetric part- 
ners). Fig. shows the dispersion relation u(k) satisfied by a fermion with thermal mass 
m at temperature T. For simplicity, we approximate it with a Lorentz-invariant relation 
uo 2 = m 2 + k 2 , shown by the solid line in fig. [1^. As discussed later, this approximation has 
a negligible impact on our results. 

Temperature corrections to the SM and MSSM particle masses are well known 
The relevant formula? are collected in appendix [B] and the numerical values of m/T are 
plotted in fig. 01 as function of the temperature (assuming a Higgs mass rrih = 115 GeV). 
The computation of thermally corrected SM decay rates and reduced cross sections is 
performed in appendix El Here we wish to discuss the most important features of the 
results. 

The processes that affect SM thermal leptogenesis are (Feynman diagrams are plotted 
in fig. EJ) 

• the decays iV — > HL and (at very high temperature when the Higgs becomes heavier 
than Ni) H — > NL (the relative reaction density is denoted as 70, see appendix [Bj) ; 

• the AL = 2 scatterings LH — > LH and LL — > HH (7^); 

• the AL = 1 scatterings involving the top quark NiL — > Q 3 U 3 ('jss) an d LQ 3 — > NiU 3 , 
LU 3 — > N\Q 3 (754) as well as their inverse reactions (which have the same reaction 
densities, up to small CP- violating corrections). We introduce AL = 1 scatterings 
involving SU(2)®U(1) gauge bosons A: N X L HA (7^), LH -> N ± A, LA -> N X H 
(lAt) and define the total scattering rates •yssj = lHs,t + lAs,t- 

Fig. shows the reaction densities 7^, •jss, 1st, ^1n, all normalized in units of Hn Nl , as 
function of the temperature and for fhi = r ■ lAm 2 ^] 1 ^ 2 = 0.06 eV and = 10 10 GeV. 

Fig. (c) show the full set of reaction densities computed including (not including) 
the effects added in this paper. In these figures we use conventions adopted in previous 

2 We do not consider the case of quasi-degenerate right-handed neutrinos. In such a case small corrections 
which break degeneracy could not be neglected. 
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Figure 4: Feynman diagrams contributing to SM thermal leptogenesis. 

papers, and plot jd, 1h s , lm, Ja s , lAt normalized in units of Hun 1 and the 'subtracted 
AL = 2 scattering rate' (see appendix EJ 7™ b normalized in units of ifn 7 . 

Decays 

The modification in 7^ is probably the most apparent feature of a comparison between 
fig. Eb and Efc, and it occurs because at sufficiently high temperature, the Higgs be- 
comes heavier than N\ and the decay N\ — > HL becomes kinematically forbidden. For 
temperatures in the range where run — tul < TtjVi < + Wl, there are no two-body 
decays involving N\ at all. At higher temperatures the Higgs becomes heavy enough for 
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Figure 5: The SM reaction densities forrhi = (Y u Yj)nv 2 /m^ 1 = 0.06 eV and = 
10 10 GeV. Blue line: decays (jd)- Red long-dashed lines: AL = 1 scatterings ('yss^ — 
jHs,t + lAs,t)- Green dashed lines: AL = 2 scatterings (jn)- 

the H — > N\L decays to be allowed, and the heavy neutrinos are produced in the pro- 
cess H «-> N\L rather than in N\ <-» HL. Including thermal masses we get jd oc T 4 at 
T ^> niN 1 - Neglecting thermal masses gave a much smaller decay rate, 7d oc T 2 , so that 
higher order AL = 1 scatterings ^fHs,t were dominant. 

AL = 1 scatterings 

A comparison of fig. Eb,c reveals other important numerical differences. There is a signif- 
icant reduction of AL = 1 scattering rates (red dashed curves), mainly 7# t , over the full 
temperature range. This comes from two different effects: 1) the top Yukawa coupling at 
high temperatures is smaller than at the electroweak scale, e.g. yt(10 10 GeV)/yt(mt) ~ 0.6; 
2) Higgs boson exchange in the t channel mediates a long-range force, giving rise to cross 
section enhanced by In m^f ttlh- This enhancement disappears when the thermal Higgs 
mass, rriH ~ 0.4 T, is used in place of the zero-temperature Higgs mass, ran ~ 100 GeV. 
Such reduction of the AL = 1 scattering rates turns out to be the most important modifica- 
tion and affects leptogenesis in two different ways: less washout of the leptonic asymmetry 
at T -C m>N 1 and slower thermalization of N± at T ^> m^. It is partially compensated by 
the inclusion of AL = 1 scatterings involving gauge bosons. 

AL = 2 scatterings 

There is a new resonant enhancement of AL = 2 scatterings LH <-» LH mediated by A^. 3 
At low temperatures the virtual N\ can be on-shell when exchanged in the s-channel, so 
that 7tv is enhanced by the s-channel resonance at s — m 2 Ni . At high temperatures the 
virtual N\ can be on-shell when exchanged in the u-channel, so that 7at is enhanced by the 

3 Scattcrings mediated by the heavier ^2,3 are included as described in section 0] 
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w-channel resonance at u — rn? Ni . For intermediate temperatures there is no resonance. 

The s-channel resonance is regulated, as usual, by the N\ width. On the contrary the 
new u-channel resonance occurs when N\ is stable, because its decays are forbidden by 
thermal masses. The divergence in the cross section is eliminated by the presence of an 
imaginary part in the N\ propagator, corresponding to the thermal damping rate caused 
by the interactions with the plasma. Although this effect depends in a complicated way 
on N\ thermal motion with respect to the plasma, in narrow width approximation (i.e. for 
small N\ Yukawa couplings) resonant enhancements give 77V — 7d/4 (see appendix lB|) . 

Therefore we performed a precise computation of the decay rate (including also Pauli- 
blocking and stimulated emission factors, as discussed in appendix |Bj) and computed 7^ b = 
In — 7d/4, the contribution to the AL = 2 scattering rate due only to off-shell scatterings. 
Indeed, this quantity enters the Boltzmann equations, because the contribution of on- 
shell Ni exchange is already taken into account by successive decays and inverse decay 
processes, LH <-> Ni «-> LH, and has to be subtracted in order to avoid double counting. 
In appendix[X]we show how to properly perform the subtractions to the Ni propagator. Our 
result differs from the one of ref . [SJ HO] • Indeed, the subtraction method used in ref . [SJ ED] 
leaves a spurious contribution which effectively double counts the decay process. This is why 
our 7™ b in fig. [SJd no longer has the "off-shell resonance" found by previous computations, 
shown in fig. Et. At leading order in the neutrino Yukawa coupling Y v iDilSsilst (X - 
(Y v Yj)u and 7^ b oc {Y V Y^)\ X . Therefore, we find that the off-shell contribution is relevant 
only when Y v ~ 1, which is not the case in our example of fig. El When Y v ~ 1 a fully 
precise computation should include also AL = LL — > N1N1 scatterings, which would play 
a minor role, affecting the Ni abundancy. More importantly, in such a case off-shell AL = 2 
scatterings suppress n# exponentially, because (unlike AL = 2 mediated by on-shell N±) 
at T <rriN 1 they are not suppressed by the Ni abundancy. 

Furthermore, fig. 03 allows to get the rates for other values of fh\ by applying the 
appropriate rescaling. For other values of m Nl the rescaling of the Yukawa couplings 
(Y u Yj) u = mim Nl /v 2 is the dominant effect, but is not the only one. One needs to 
recompute the rates taking into account the running of the couplings. 

3.3 CP violation at finite temperature 

In this section we investigate the finite-temperature effects on the CP asymmetries 

£ = 7 cq (^/)-7 cq (^7) f28l 
1 7 eq (^/)+7 eq (^7) ' 

where 7 cq are the thermally averaged decay rates. The decay processes relevant to our 
analysis are N\ — > LH, which is allowed for > m^(T) + m^(T) and H — > LN±, which 
is allowed at higher T, when mjy(T) > m^iT) + and is CP- violating only because of 
purely finite-T effects. 

The issue of CP- violating decays at finite T was already investigated in ref. [2E], although 
neglecting thermal masses. As we will see, the effect of the masses is crucial, giving a 
non-trivial T dependence of the CP violation. The effect of thermal masses is taken into 
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account by using the one-loop finite temperature resummed propagators, and by using 
modified dispersion relations, as discussed in section |21 

We choose to work in the rest frame of the plasma, where = (1, 0, 0, 0). In this way 
the finite-T Feynman rules are simplified, while the kinematics is more complicated, since 
we have to consider decaying particles in motion. Although we have performed our calcu- 
lation in the general case, neglecting the thermal motion of the decaying particles allows 
to write reasonably accurate analytical approximations. As already said, we approximate 
the complicated dispersion relation for fermions with uj 2 = k 2 + m 2 F . 

The CP asymmetries in the relevant decays come from interference between the tree- 
level decay amplitude with the one-loop contributions. There are two relevant one-loop 
diagrams: the so-called vertex and wave- function contributions (shown in fig. ITd^I . We 
compute the imaginary part of the one-loop graphs using the Cutkosky cutting rules at 
finite temperature [27], which are more complicated than at T = (even in the absence 
of type-2 vertices). While at T = most cuttings (in our case two of the three possible 
cuttings) give no contribution due to energy-conservation, this is no longer true at T ^ 0, 
since particles may absorb energy from the plasma. Formally, this means that we must 
also consider cut lines with negative-energy on-shell particles. 

Nevertheless, we ignore cuts which involve heavy right-handed neutrinos A^,3, because 
they are suppressed by a Boltzmann factor exp(— 3 /T), which is negligibly small since 
we assume a hierarchical spectrum m^ 2S 3> m^, and we work at T much smaller than 
mN 23 - Therefore, as illustrated in fig. EH we can restrict ourselves to the standard cutting 
of the Higgs and lepton lines, but with energy flows in both directions. 

CP-asymmetry in Ni decay 

We first consider the CP asymmetry in N± decay, 



The full result employed in our leptogenesis code is presented in appendix Here we 
present a simple analytic approximation obtained neglecting the thermal motion of Nx 
with respect to the plasma, which is justified at T < and still reasonably accurate at 
higher temperatures (see fig. EK) • The result is 



7 cq (A^! -> HL) - 7 eq (iVi -> HL) 
7 eq(JVi ->• HL) + 7 ec i(iV 1 -> HL) ' 



(29) 





(30) 



The function / describes the usual result at T = which, in the SM, is given by 



f(x) = Vx 



x 



x-2 



- 1 



(1 + x)ln 



( 



1 + x 



x 



) 



3 



(31) 
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The thermal correction is described by the function R t , given by 4 

(32) 



Re = 16 [u(l + a L ) + b L ]J[l + f H -h- 2f H f L ] J 

m Nl 



dS H /d(jj dS L /duj 
dS H /dk d5 L /dk 



where ||M|| = | det M\, Kl = (uj, k) is the energy-momentum quadri- vector of L, and the 
functions = a{Ki) and = b(Ki) are defined in eq.s ()16I17|) . The Fermi-Dirac and 
Bose-Einstein distributions 

h = (^ /T + 1)" 1 f H = {e E »' T -i)- x 

in the third factor are evaluated at the fermion energy El = uj and at the boson energy 
Eh = jvj — uj. The first term is obtained by computing the relevant Feynman graph in 
the limit m^ 2Z 3> m^ 1 and dividing by the tree level rate. The Jacobian J is obtained 
when imposing the on-shell conditions for the cut particles H and L 

5 H = (m Nl -Lu) 2 -k 2 -m H = 5 L = [(1 + a L )uj + b] 2 - (1 + a L fk 2 = (33) 

which fixes the values of u and k in terms oim Nl , m L , m H . As discussed in section El the 
equation for L has two different solutions: 'particles' and 'holes'. A numerical computation 
shows that the 'hole' contribution to the CP-asymmetry is negligible because, as explained 
in section |21 relativistic holes have negligible interactions. The 'particle' contribution is 
well approximated by inserting in eq. (}3*2|) the values of uj and k 



m 2 N + mi — mi 



U = 2^ ' k = - < ^ 

obtained approximating 5l ~ uj 2 — k 2 — m\ = 0. 

The main feature shown in fig. El is that e(T) goes to zero as the temperature increases 
and the process becomes kinematically forbidden. This happens because the particles in 
the final state coincide with the cut particles in the loop: L and H. Therefore the threshold 
at which the cut particles can no longer be on the mass-shell is the same at which the decay 
becomes kinematically forbidden, i.e. when ~ + m^. 

There is another important effect which gives an additional suppression. The 1 + fs — 
fp — 2f B f F factor in eq. (J3*2~|) was first derived by the authors of ref. who, neglecting the 
L and H thermal masses and thus setting uj = found it to be equal to 1. However, 

only if the arguments of the Bose-Einstein and Fermi-Dirac distributions are the same, 
there is a peculiar cancellation: ]b — ]f — 2]b]f — 0. Physically this cancellation can be 
understood as a compensation between stimulated emission and Pauli blocking. Only if 
bosons and fermions enter with the same energy, an exact cancellation holds. 



CP-asymmetry in H decay 

The computation of the CP asymmetry in Higgs decay, 



= 7 eq ( ff^ iVl L)- 7 ^(g^iV 1 £) 



4 If -ZV2.3 are not much heavier than N± one can easily include effects suppressed by higher powers of 
ti JVi /tti jv 2 3 • obtaining a more lengthy analytical expression. 
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Figure 6: Thermal corrections to the CP asymmetry. e(T)/e(T = 0) as a function 
of temperature for m?q 2 3 3> = 10 10 GeV in the SM (left plots) and in the MSSM (right 
plot, the dots indicate the various thresholds for Ni and N\ decays). In the SM plot the solid 
line shows our more accurate result derived in appendixlTX while the dashed line shows the 
approximate result described in the main text, obtained neglecting the Ni thermal motion. 



is similar to the previous one, although with some important differences. Also in this case 
there are two relevant cuts (in wave and vertex one loop diagrams, see fig. 1X7)1 which involve 
the Higgs and lepton lines. The difference is that such graphs would have no imaginary 
parts with the usual Feynman rules at T = 0, as cuttings of H and L would be kinematically 
forbidden. On the contrary, at finite T absorption of particles by the plasma allows also 
negative energies in the cuts. The asymmetry e# turns out to be proportional to the 
purely thermal factor fjj — f^ — Ifufh- The non-standard cut (see fig. 117)) implies a more 
complicated kinematics which does not allow us to obtain an analytic result for e#. The 
computation is presented in appendix where we neglect effects due to particle motion 
and due to non-trivial fermion dispersion relation, since eventually e# turns out to have 
a negligible effect on the final results for leptogenesis. e# approaches a constant value at 
high T 3> m Nl (see fig. EJ). 



CP-asymmetries in the MSSM 

The situation in the MSSM becomes more complicated than in the SM because we must now 
study both N\ and N\ decays, each having two possible decay channels, with each channel 
having more diagrams. The CP-asymmetries at zero temperature were first computed 
correctly in ref. [21] • We do not study H and H decays, which appear only at high 
temperature. 

The computations at finite temperature are analogous to the SM case and described 
in appendix [Dj The result is different because the cut particles in the loop are different 
from the final-state particles: the threshold at which the cut particles can no more go 
on-shell is different from the one at which the decay becomes forbidden. Therefore, in the 
single decay modes, CP-violation disappears either before or after the decay mode becomes 
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kinematically forbidden. 

The source of CP-asymmetry which appears in the MSSM Boltzmann equations are 
the CP-asymmetries ejVi and in Ai and Aq decays, averaged over the different decay 
channels Aq — > HL,HL and Aq — ► HL,HL. As illustrated in fig. Eb £n 1 and behave 
in rather different ways. 

Like in the SM, goes to zero at the same temperature threshold at which Aq decays 
become kinematically forbidden, as a consequence of the fact that both decay channels 
contribute to the imaginary part of both decay modes. 

On the contrary, for , only one-loop diagrams with internal bosons H, L contribute 
to the CP-asymmetry of decays into fermions H, L and viceversa: since thermal corrections 
make bosons heavier than the corresponding fermions (see appendix El) vanishes when 
N\ decays into fermions are still kinematically allowed. The dashed line in fig. Efc> shows 
the final result. Thermal corrections are significant at T ~ m^, but give almost no effect 
at lower temperatures. This happens in a non-trivial way. Decays into scalars have a rate 
significantly enhanced by stimulated emission and a CP-asymmetry significantly suppressed 
by Pauli blocking, while the opposite happens for decays into fermions. If thermal masses 
can be neglected, the two effects compensate each other, as noticed in ref. [2Ej- This 
cancellation no longer takes place when thermal masses become sizable, giving rise to the 
behavior of (T) shown in fig. Eb- 



4 Leptogenesis in the Standard Model 

We assume that right-handed neutrinos are hierarchical, rriN 2a 3> so that we have to 
study the evolution of the number density of Aq only. In such a case the final amount of 
B — C asymmetry Yb-c = n B-cl ' s generated by Aq assuming no pre-existing asymmetry 
can be conveniently parameterized as 

Yb-c = -ejvi V Y^(T > m Nl ). (36) 

Here is the CP-asymmetry parameter in Aq decays at zero temperature, and Y^(T ^> 
mjvi) = 135C(3)/(47r 4 g*), where g* counts the effective number of spin-degrees of freedom 
in thermal equilibrium (g* = 106.75 in the SM with no right-handed neutrinos) 5 , r] is an 
efficiency factor that measures the number density of Aq with respect to the equilibrium 
value, the out-of-equilibrium condition at decay, and the thermal corrections to e^. Re- 
calling that, after reprocessing by sphaleron transitions, the baryon asymmetry is related 
to the B — C asymmetry by 

s_ ~ 66 + \?>n H s ' [ ' 

J Thc formula used in our numerical code includes leading order thermal effects from quarks, leptons 
and gauge bosons (neglecting Yukawa couplings): 



p = pr + Pm, pr 



427 n 2 7 2 19 2 25 2 
~30 ~ 4 53 ~ 32 52 ~ 96 5y 



T 4 a - APR 
3T ' 
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where is the number of Higgs doublets, for the SM we find 

— = -1.38 x 10~ 3 e Nl rj. (38) 

s 

Assuming the 'standard' ACDM cosmological model, BBN and WMAP measurements 
imply 

— = (6.15 ± 0.25) x 10" 10 with s = 7.04 n T (39) 

Tiny 

Computing rj is the most difficult part of the calculation, since it is obtained from 
numerical solution of Boltzmann equations. In general the result depends on how the 
lepton asymmetry is distributed in the three lepton flavours. For simplicity one usually 
ignores flavour issues and solves the approximated Boltzmann equation for the total lepton 
asymmetry described in appendix El If all mixing angles of left and right-handed neutrinos 
are large, as data might suggest, this 'one-flavour approximation' is accurate up to 0(1) 
corrections. In order to include flavour factors one must solve the Boltzmann equations for 
the 3x3 density matrix of lepton doublets, as discussed in ref . j2H] • 

The AL = 2 scatterings mediated by heavier right-handed neutrinos m^ 2 3 3> m^ 
are relevant at > 10 14 GeV. Below, they can have 0(1) effects if neutrinos are quasi- 
degenerate. Although these scatterings are produced by the same effective dimension 5 
operator that generates neutrino masses, their contribution to &Ns,t depends on unknown 
high-energy parameters: the flavour composition of the neutrino coupled to Aq. There- 
fore, following ref. [2E], we introduce a parameter £ that, in eqs. E2HHS1 parameterizes the 
unknown contribution of A^, which is important only if > 10 14 GeV. If there were 
only one neutrino flavour with mass m u , the value of £ would be £ = m v /rh\ — 1. With 
three neutrinos we cannot even write the precise definition of £, as flavour factors cannot 
be correctly included in Boltzmann equations valid in 'one-flavour approximation'. In our 
numerical results we assumed £ = max(l, m atm /mi): if all mixing angles are large this 
choice is reasonably correct in all the parameter space. 

Having fixed £, in 'one-flavour approximation' rj depends only on two parameters: 

mi = (Y v Yj)iiV 2 /m Nl and m Nl 

(and almost only on mi if <C 10 14 GeV i.e. when AL = 2 scatterings mediated by 
Aj off-resonance are negligible). This well-known fact remains true also when thermal 
corrections are included, as can be seen from appendix [Bl The parameter mi is 'the 
contribution to the neutrino mass mediated by Aq'. To see what this means in practice, 
let us temporarily assume that the left-handed neutrinos have a hierarchical spectrum 
mi <C m 2 = m sun <C m 3 = m atm = jAm^] 1 / 2 . In such a case mi = m atm if Aq exchange 
gives rise to the atmospheric mass splitting, or mi>m sun if Aq gives rise to the solar 
mass splitting. A smaller mi > mi can be obtained if Aq gives rise to m 1; which can be 
arbitrarily small. A mi larger than m atm can be obtained if A^ exchange cancels out the 
Aq contribution to neutrino masses. Stronger (weaker) restrictions are obtained if there 
are less (more) than 3 right-handed neutrinos. 

In conclusion, measuring neutrino masses does not fix mi and m^, which remain as 
free parameters. Therefore we compute r] as function of and of mi renormalized at 
the high scale m^ (at high scales mi is about (20 -r 30)% larger than at low energy). 
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Figure 7: Evolution of (blue curves) and \Yb-lI 'ejvil ( re d curves) with temperature 
in the SM. We fix mjvj = 10 10 GeV and miim^) = 0.06 eV. T/ie dashed line shows the 
thermal abundance ofY^. Left plot: full computation, the efficiency is 77 = 0.0036. Right 
plot: no new effect included (and on-shell scatterings incorrectly subtracted), t] = 0.0017. 



Results 

Figure [3 shows the evolution of the N\ and B — £ abundances at our sample 'atmospheric' 
point: rhiirriN^ = r(Am^ tm ) 1 / 2 = 0.06eV and mjvi = 10 10 GeV. For these values the N\ 
abundancy remains close to thermal equilibrium, so that leptogenesis is mainly determined 
only by the later stages of the evolution at relatively small temperatures. This explains 
why, despite the significant variations at higher temperature, there is only a mild correction 
to the final baryon asymmetry. Proper subtraction of on-shell scatterings reduces wash- 
out by a 3/2 factor. This gives a 3/2 increase of the efficiency, as can be seen from the 
analytical approximation of ref. |29j . 

We now present our results for the thermal leptogenesis efficiency parameter 77. We 
assume no pre-existing B — C asymmetry, and we study the three cases of 

(0) zero initial iV\ population, Y^ l = at T > mjy r This case can be realized e.g. if an 
inflaton field reheated the universe decaying mostly into SM particles. 

(1) thermal initial N\ population, Y^ 1 = Y^ at T ^> mjvi- This case can be realized in 
presence of extra interactions at T ^ m,N 1 , mediated e.g. by a heavy Z' boson related 
to SO(10) unification. 

(00) dominant initial Ni abundancy, Pn x ^ Pr at early times. This case can be realized 
e.g. if an inflaton field reheated the universe decaying mostly into Ni. 

In order to study the latter case we modified the Boltzmann equations employed in previous 
analyses including the effects of Ni reheating (in the case (1) Ni reheating is relevant only 
if m x < 10" 6 eV). 

If m^Vj -C 10 14 GeV the efficiency parameter 77 depends almost only on fh\. In fig. |H]we 
show 1 7] I as function of fhi for = 10 10 GeV. If fh\ > 10 -2 eV neutrino Yukawa couplings 
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Figure 8: Efficiency n of leptogenesis in the SM, assuming zero (dashed red line), 
thermal (continuous blue line) or dominant (long dashed green line) initial Ni abundancy. 
Upper plots: n as function of m\ (renormalized at m^) for = 10 10 GeV. Lower plots: 
contours ofr)(rhi, m^) = lO -6 ' -5 ''"' 0,1 . In the shaded regions the neutrino Yukawa couplings 
are non-perturbative. 

keep N\ so close to thermal equilibrium that 77 does not depend on the unknown initial N\ 
abundancy (similarly, an eventual pre-existing lepton asymmetry would be washed out if 
N\ Yukawa couplings act on all flavours). 

At smaller ih\ the efficiency 77 depends on the initial Ni abundancy, ranging between 
the limiting cases (0) and (00), as illustrated by the gray band in fig. |HJ As expected, the 
maximal value of n ~ g* is reached at fhi ~ fh\ = 256 v /^f 2 /3Mpi = 2.2 x 10~ 3 eV in 
case (00). In such a case, rj decreases at « m], because N\ decays out-of-equilibrium 
at temperature ~ m^^fh\jm\ €%! so that N\ reheating washes out some lepton 
asymmetry. In more physical terms, the particles H, L emitted in N\ decays have energy 
larger than the temperature T, and split up in ~ m Nl /T^ { particles without correspond- 
ingly increasing the lepton asymmetry, so that n ~ g* ^m\/rhi. 

When 771 1 < 10~ 6 eV, N\ reheating starts to be significant even in case (1) giving 77 < 1. 
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Figure 9: Allowed range of fh\ and mjvj /or leptogenesis in the SM and MSSM assuming 
r«3 = max(mi,m atm ) and £ = mz/fhi. Successful leptogenesis is possible in the area inside 
the curves (more likely around the border). 

In fact, even if Ni initially has a thermal abundancy Pnx/ Pr ~ 9 Nil 1 9* "C 1, its contribution 
to the total density of the universe becomes no longer negligible, PnJ Pr ~ (gN x m N^)l (9*T), 
if it decays strongly out of equilibrium at T C m^. For the reasons explained above, this 
effect gives a suppression of rj (rather than an enhancement), and for very small m\ the 
case (1) and (oo) give the same result. 

The lower panel of fig. |H1 contains our result for the efficiency \rj\ of thermal leptogenesis 
computed in cases (0), (1) and (oo) as function of both fh\ and mjvi- At mjvj ^ 10 14 GeV 
non-resonant AL = 2 scatterings enter in thermal equilibrium strongly suppressing rj. 
Details depend on unknown flavour factors. 

Our results in fig. |H] can be summarized with simple analytical fits 



3.3 x 10" 3 eV 



mi 



+ 



mi 



0.55 x 10" 3 eV 



1.16 



in case (0) 



(40) 



valid for 10 14 GeV. This enables the reader to study leptogenesis in neutrino mass 

models without setting up and solving the complicated Boltzmann equations. 



Implications 

Experiments have not yet determined the mass m 3 of the heaviest mainly left-handed 
neutrino. We assume m 3 = max(m 1 , m atm ). Slightly different plausible assumptions are 
possible when fhi m m atm , and very different fine-tuned assumptions are always possible. 
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Figure 10: Leptogenesis bound on neutrino masses. The plot shows the measured 
baryon asymmetry (horizontal line ) compared with the maximal leptogenesis value as func- 
tion of the heaviest neutrino mass m 3; renormalized at low energy. Error bars are at 3a. 

As discussed above, we assume £ = m^/mi (the parameter £ controls AL = 2 scatterings 
mediated by A^). These assumptions do not affect the absolute bounds on the masses of 
left-handed and right-handed neutrinos that we now discuss, but allow to present them in 
one simple plot, fig. El 

The crucial assumption behind fig. El is that right-handed neutrinos are very hierarchical. 
Under this hypothesis the CP asymmetry is bounded by the expression given in [22] (see 
also (201 EH), that in the hierarchical and quasi-degenerate light neutrino limits simplifies 
to 

3 rriN 1 (m 3 — mi) ( \ — m\jm\ if m\ <C friz 



\e Nl \<T2 1 5 x \ F\ 27-^2 -t ■ ( 41 ) 

167T v 2 y\ — m\jm\ if mi ~ 777.3 

where all parameters are renormalized at the high-energy scale ~ m^ 1 and m\ = m\ + 
Am^ tm + Amg Un . The 3cr ranges of m atm and of n B /n 7 imply the lower bound 

45xl0 8 GeV f 2.4xl0 9 GeV in case (0) 
m Nl > — — ^JfZl. > ) 4.9 x 10 8 GeV in case (1) (42) 
71 { 1.7 x 10 7 GeV in case (00) 

with the absolute bound realized in case (00). 

The allowed regions are shown in fig. El as function of rh\. The bound on mj$ x becomes 
stronger if left-handed neutrinos are heavier than what suggested by oscillations, m 3 > 
m atm , until thermal leptogenesis can no longer generate the observed ng (see fig. EJ) giving 
an upper bound on the mass of degenerate neutrinos (renormalized at some unspecified 
scale j2Hj) of about 0.1 eV [Toj . This happens because 77 decreases with rhi in the region of 
interest and because the bound in eq. (}4*l"l) on becomes stronger when neutrinos become 
heavier. 6 



6 The bound on neutrino masses is saturated around ~ 10 13 GeV and fh\ ~ 0.1 eV. For these values 
computing leptogenesis in 'one flavour approximation' is reliable |22| . We recall that we are assuming that 
the Higgs quartic coupling remains relatively small up to high energies. 
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In order to study precisely the bound on neutrino masses we relax our simplifying 
assumptions on £ and m 3 and compute the maximal baryon asymmetry generated by 
thermal leptogenesis as function of m^,. The results is shown in fig. including the new 
effects discussed in this paper and combining errors on Am^ tm and on the baryon asymmetry 
in quadrature we get m u < 0.15 eV at 99.73% CL (i.e. 3a). The small difference with respect 
to previous results is due to various factors: correct subtraction of on-shell scatterings 
(makes leptogenesis 50% more efficient and the bound on neutrino masses 7% weaker); 
renormalization of neutrino masses (makes the bound 7% stronger); renormalization of X t , 
inclusion of gauge scatterings, of thermal corrections, updated experimental determination 
of Am^ tm and of rig/n^, revised (weaker) upper bound on the CP-asymmetry [2*2*] . 

We stress that the bound in eq. (j41j) on the CP-asymmetry holds because we assumed 
mjvi C ffljv 23 : if right-handed neutrinos were instead quasi-degenerate CP violation in 
mixing can give an arbitrarily large CP-asymmetry, e^ x ~ 1 and all bounds that we dis- 
cussed evaporate. In particular the leptogenesis bound on the neutrino masses holds under 
the dubious assumption that hierarchical right-handed neutrinos give quasi-degenerate left- 
handed neutrinos. 

4.1 Simple approximation 

While in our calculations we use our full code, we now also discuss how the computation 
can be significantly simplified by including only the following new ingredients which, a pos- 
teriori, turn out to be numerically most relevant. Neglecting all these effects we reproduce 
the results in ref. [HI]. Effects 2, 3, 4 were included in ref. j2H] : neglecting the other ones 
we reproduce their results. Effect 2 was mentioned in ref. |2*3*] . 

1 Proper subtraction of the on-shell JVi propagators in AL = 2 scattering processes, as 
discussed in Appendix A. 

2 Neutrino masses have to be renormalized at the proper energy scale ~ m^ 11 both 
when computing r\ and e^. 

Typically these effects give 0{\) corrections. For fhi > 10~ 3 eV the efficiency increases by 
up to almost a factor 2, due to the suppression of AL = 1 wash-out processes caused by 
the following thermal effects: 

3 Temperature corrections to the Higgs boson mass, m# ~ 0.4 T, must be included at 
least when computing the IR-enhanced LN — > Q3U3 interaction rate. 

4 The top Yukawa coupling must be renormalized at the proper energy scale ~ m^. 
These variations are partially compensated by 

5 Inclusion of previously neglected scatterings involving weak gauge bosons A, N\L <-» 
HA, N\H <-» LA and N\A «-> HL. These extra scatterings are sizable because 
#2 > A 4 at energies above 10 9 GeV: \ t is no longer the dominant coupling. These 
processes have been recently considered for the time in ref. jH2] 
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Figure 11: Efficiency rj of leptogenesis in the MSSM, assuming zero, thermal or dom- 
inant initial A r 1 ,A r 1 abundancy. In the shaded regions the neutrino Yukawa couplings are 
non-perturbative. 

Including these contributions one gets an excellent approximation for rh\ > 10~ 3 eV. In 
order to get an approximation which is accurate also at rh\ <C 10~ 3 eV, in case (0) one 
needs to include one more effect: 

6 Thermal corrections to the CP- violating parameter e^. 

which turns out to have a sizable impact due to a more subtle reason. 

Neglecting washout scatterings (which are small at rhi <C 10~ 3 eV) and the temperature 
dependence of the CP asymmetry, in case (0) the Boltzmann equations are solved by 
Yg-£ = +ejViYivi- In this approximation the lepton asymmetry generated in inverse-decay 
processes at T ^> m Nl when Y Nl < is exactly cancelled by the asymmetry generated 
later in N\ decays when Y^ 1 > Y^. Consequently the effect which dominantly breaks 
this cancellation has a numerically important impact even if it is 'small'. 1) Washout 
interactions, dominated by 7/5, erase the lepton asymmetry generated at earlier stages more 
than the one generated at later stages, giving a small positive r]. 2) Thermal corrections 
to e slightly reduce it at small temperatures T ~ m^/A, giving a small positive rj. 7 

Finally, if one wants to study cases where N\ can give a substantial contribution to 
the total energy density, one must include this correction in the Boltzmann equations, as 
described in appendix 1X1 

7 The solid line in fig. [5^ shows our most accurate result for thermal corrections to e, that we employ 
in numerical computations. The enhancement at T ~ 0Arri]\i 1 comes from the quantum statistics factor, 
1 + fn — fi, — that can be larger than one when both thermal masses and N\ motion with respect 

to the plasma are taken into account. In case (0) for vn\ 10 -3 eV the Ni energy spectrum deviates from 
the thermal Fermi-Dirac distribution that we assumed. Using the slightly less accurate thermal correction 
to e, obtained neglecting thermal motion of iVi (dashed line in fig. EJi) typically affects the final result by a 
O(l) factor. Other minor effects might be important: the exact dispersion relation at finite temperature, 
thermal corrections to couplings, higher order corrections, and CP-violation in AL = 1 scatterings. 



23 



5 Leptogenesis in the MSSM 



In the case of the supersymmetric extension of the SM, the computation becomes more 
involved because of the presence of many new particle degrees of freedom. In particular, 
a lepton asymmetry is generated in decays of both the right handed neutrino N% and the 
right-handed sneutrino N\. Since supersymmetry breaking can be ignored (except in spe- 
cial cases |TT) I12|). the computation is somewhat simplified by the identities demanded by 
supersymmetry, such as = m N 1 and = . Thermal corrections break supersym- 
metry, so that it is rather ponderous to include them. At the moment, we do not attempt 
to make a full calculation in the supersymmetric case. 

On the other side, a full computation might be not necessary. Including only the thermal 
effects which in the SM turn out to be dominant as discussed in section l4~T| could be a good 
approximation also for the MSSM. Therefore we do not compute all relevant cross sections 
including finite temperature effects, but we adopt the MSSM cross sections of ref. |33j . 
inserting the temperature-dependent top Yukawa coupling and Higgs boson mass in IR- 
enhanced processes, and performing a correct subtraction of on-shell resonances (which 
again reduces ^-mediated washout scatterings by a 3/2 factor). Thermal corrections to 
N\ and N\ decays and their CP-asymmetries are computed in appendix O neglecting the 
thermal motion of Ni, N\ (in the SM case this would not be a very good approximation 
at small mi). Finally, we do not include AL = 1 scatterings involving gauge bosons and 
gauginos. 

Our MSSM results have been obtained under these approximations, assuming moder- 
ately large values of tan/5 ~ 10. Low-energy thresholds make the top Yukawa coupling at 
high energy uncertain by about a factor of 2, as discussed in ref. |24j . 

Proceeding as in the case of the SM, we find that the asymmetry generated in the 
MSSM is 

Y B „c = -i? e Nl + 1£) (T » m Nl ) , (43) 

where is the neutrino (or sneutrino) CP-asymmetry at low temperature (equal for lepton 
and slepton final states) [2 



g{x) = -y/X 



±Y^m g (<\ (44) 



2 , fl + x 
+ In 



x — 1 \ x 



(45) 



x 



The number of effective degrees of freedom in the MSSM without right-handed neutrinos 
is g* = 228.75. Using eq. (|3Tj) with n H = 2, we obtain 



s 



-1.48 x lO^e^. (46) 



The MSSM results, analogous to those obtained from the SM and previously discussed, 
shown in figs. IUUTT1 are similar to their corresponding SM results. It is difficult to compare 
with previous results, which have not been presented in a systematic way. 
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zero Nj 



thermal Nj 




Figure 12: Regions o/fm^mjvj) plane where soft leptogenesis can produce the observed 
baryon asymmetry for ImA < TeV and ^Bm^ 1 = 100 GeV (solid line), 1 TeV (dashed), 
lOTeV (long-dashed) We have assumed a vanishing (left) or thermal (right) initial sneu- 
trino density. 



Extra signals may come from lepton flavour violating decays like /i — > cy, induced in 
supersymmetric see-saw models by the (unknown) neutrino Yukawa couplings [HI]. Some 
predictive minimal models allow to predict these rates in terms of the measured neutrino 
masses and baryon asymmetry. For example, using our revised leptogenesis computation, 
the prediction of ref. jSH] for BR(/i — > ej) gets lowered by one order of magnitude. 



Soft leptogenesis 

"Soft leptogenesis" [TT) IT2*| is a supersymmetric scenario of leptogenesis which requires 
only one heavy right-handed neutrino. The interference between the CP-odd and CP-even 
states of the heavy scalar neutrino resembles very much the neutral kaon system. The mass 
splitting as well as the required CP violation in the heavy sneutrino system comes from 
the soft supersymmetry breaking A and B terms, respectively associated with the Yukawa 
coupling and mass term of N\. The non- vanishing value of the generated lepton asymmetry 
is a pure thermal effect, since at T = the generated lepton asymmetry in leptons exactly 
cancels the one in sleptons. 

Here we improve the results of ref. ^2] , taking into account thermal, Pauli blocking and 
stimulated emission corrections to N%,N\ branching ratios (which significantly enhance 
decays into bosons) to their CP asymmetries (computed in appendix O including thermal 
corrections, but neglecting the thermal motion of Ni,Ni). We recall that in this scenario 
e N->LH nas an opposite sign to and 



e ^ te( T = o) = - e ^(r = o) = i ^-^ (47, 
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We refer to ^2] for further details, and present our improved results in fig.ll2[ where we plot 
the regions in (m^mjvj plane where successful leptogenesis is possible. We have assumed 
vanishing (left) or thermal (right) initial sneutrino density, ImA < TeV and ^BrriN 1 = 
100 GeV (solid line) and 1 TeV (dashed line), 10 TeV (long dashed line). 'Soft leptogenesis' 
needs a S-term much smaller than what suggested by models of supersymmetry breaking, 
but not unnaturally small ^2]- Soft leptogenesis allows to lower the bounds on ttin 1 - 

6 Leptogenesis with reheating of the universe 

By now there is a wide consensus that the early universe underwent a primordial stage of 
inflation [HE] responsible for the observed homogeneity and isotropy of the present universe, 
as well as for the generation of the cosmological perturbations. 

The radiation-dominated era of the universe is usually assumed to be originated by the 
decay of the coherent oscillations of a scalar field, the inflaton field, whose vacuum energy 
has driven inflation For such a reason the reheating process is often associated with the 
final stage of inflation. However, we point out that reheating could have been episodic, 
with several reheat events after inflation. We will be interested in the final reheating which 
may just as well have been the result of the decay of a weakly coupled scalar field unrelated 
to inflation, for instance a modulus. For this reason the scalar field 0, whose decay leads 
to reheating, will not be referred to as the inflaton. 

The decay of the scalar field <p into light degrees of freedom and their subsequent 
thermalization, called reheating, leaves the universe at a temperature Trh, which represents 
the largest temperature of the plasma during the subsequent radiation-dominated epoch, 
when temperature is a decreasing function of time. The onset of the radiation dominated 
era is in fact placed at the temperature Trh, i-e. at the end of the reheating phase. 

Usually Trh is assumed to be very large and this is the assumption we have made in the 
previous sections. However the only information we have on the smallest value of Trh is 
from requiring a successful period of primordial nucleosynthesis, Trh ^ 1 MeV. Therefore, 
from a phenomeno logical point of view, Trh is actually a free parameter 8 . Any scenario of 
baryogenesis based on the out-of-equilibrium decay of some heavy particle depends crucially 
on the assumption that these particles were generated during the reheating process with 
abundances sufficiently large to generate the observed baryon asymmetry. During reheat- 
ing, particles are generated through thermal scatterings and quickly thermalize. Among 
them, right-handed neutrinos may be also produced but their number depends strongly 
on the reheating temperature. If the latter is too small, the thermal bath does not give 
rise to a number of right-handed neutrinos large enough to produce the observed baryon 
asymmetry. This leads to a lower limit on Trh- Computing this lower bound is the goal of 
this section. 

During the reheating epoch, the energy density of the universe is dominated by the 
coherent oscillations of a scalar field <fi. Considering for the moment the case of small 

8 Low reheating scenarios lead as well to a new perspective on baryogenesis to the possibility that 
massive neutrinos may play the role of warm dark matter [SHI or to a change in the predictions of the relic 
abundance and resulting model constraints of supcrsymmctric dark matter, axions, massive neutrinos, and 
other dark matter candidates pHU Ml ■ 
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abundance of right-handed neutrinos 9 , as a first step we assume that the dynamics of 
reheating is described by the Boltzmann equations for the energy densities p^R of the two 
coupled components: the unstable massive field (ft and the radiation R |HJ 021 EH] 



dt 
dpR 
dt 



-3Hp^ - I> , (48) 
-AH p R + , (49) 



where H = a /a = ^/87r(p0 + p^)/3M| 1; and Mpi is the Planck mass. Here we have 
assumed that the relativistic decay products of the scalar field rapidly thermalize and form 
a relativistic bath of temperature T (for a discussion about this point see Ref. [12] )• The 
key point of our considerations is that reheating is not an instantaneous process. On the 
contrary, the radiation-dominated phase follows a prolonged stage of matter domination 
during which the energy density of the universe is dominated by the coherent oscillations of 
the field 0. The oscillations start at time Hj 1 and end when the age of the universe becomes 
of order of the lifetime T^ 1 of the scalar field. At times Hj 1 < t < T^ 1 the dynamics of the 
system is quite involved. During this stage the energy density per comoving volume of the 
(f> field decreases as exp(— T^t) and the light decay products of the scalar field thermalize. 
The temperature T of this hot plasma, however, does not scale as T oc a -1 as in the 
ordinary radiation-dominated phase (where a is the Friedmann-Robertson-Walker scale 
factor) jJU 02], but reaches a maximum T M ax ~ (HiMpiY^T^ and then decreases as 
T oc a -3 / 8 , signalling the continuous release of entropy from the decays of the scalar field. 

In fact, until t <C T^ 1 assuming p^ ^> pr the system approximately evolves as 

PS) = P^(0)(a /a(t)) 3 e- r ^ (50) 

PR {t)^^ « ^r,M Plv ^ [i - (V«) 5/2 ] • (51) 

This scaling continues until the time t ~ T^ 1 , when the radiation-dominated phase starts 
with temperature T ~ T RH , defined as 



RH 



TlM? 



4ir 3 g* (T RH 



pi 



n l/4 

(52) 



The process is described by two extra parameters, ^^(0) and T^. It is convenient to replace 
them with the maximal and reheating temperatures, Tmax and Trh- Before reheating 
is completed, at a given temperature, the universe expands faster than in the radiation- 
dominated phase. Notice that Trh is not the maximum temperature during the reheat 
process. On the contrary, Tmax can be much larger than Trh- 

During reheating right-handed neutrinos may be produced in several ways. They can 
be generated directly through the scalar field perturbative decay process [IB] (this requires 
that the mass of the 0-field is larger than m^) or through nonperturbative processes taking 
place at the preheating stage j3H EH] . These mechanisms, however, introduce new unknown 

9 The case of non-negligible Ni density is discussed in appendix A. 
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parameters such as the coupling of right-handed neutrinos to <p. In this section we take a 
more conservative point of view and we limit ourselves to the case in which right-handed 
neutrinos are produced by thermal scatterings during the reheat stage, so that in the limit 
2rh ^ m 7Vi we obtain the case (0) studied in section 0] (in the opposite case of dominant 
inflaton decay into JVi, the Trh ^ Wjvi limit is given by case (oo) of section HJ. For the 
sake of concreteness we will focus on the leptogenesis scenario, but our findings can be 
easily generalized to any out-of-equilibrium scenario. Furthermore, we assume that the 
mass of the inflaton field is larger than the reheating temperature Trh; for a discussion of 
the opposite case, see ref. 



We now generalize the Boltzmann equations for thermal leptogenesis including reheat- 
ing. In the standard case it is convenient to write the Boltzmann equations that dictate the 
time evolution of the number densities rix{t) of any species X in terms of Yx(z) = nx/s 
where z = m^/T. In fact, while particle densities nx strongly depend on t because of 
the expansion of the universe, their ratios with respect to the entropy density s remain 
constant in the absence of interactions. Since the temperature T is a monotonic decreasing 
function, one usually replaces the time t with T. 

These two statements no longer hold during reheating. Nevertheless, we still find con- 
venient to write the Boltzmann equations for leptogenesis in terms of Yx(z). The first 
pre-heating phase (when T grows from zero to Tmax) is so fast that it gives no contribu- 
tion to leptogenesis: the interesting dynamics happen during the second reheating phase, 
when the temperature decreases in a non standard characteristic way, T oc a~ 3 ^ 8 . Therefore 
corrections to leptogenesis are fully described by a single parameter, the reheating temper- 
ature Trh, as long as Tmax is sufficiently larger than rajvi- 

Since the temperature is defined in terms of the radiation density by eq. (|49ft. we can 
write 

— = -AHZp R — = HZz— Z = 1 - T ^ = -—— (53) 
dt dpn dz ^Hp R T da 

Z vanishes when the maximal temperature Tmax is reached, then Z ~ 3/8 during reheat- 
ing, and finally Z ~ 1 in the standard radiation-dominated phase. Apart from this 0(1) 
correction, reheating affects leptogenesis in two main ways 1) H has a non-standard expres- 
sion: induces a faster expansion 2) <fi decays create additional matter. The Boltzmann 
equations for leptogenesis are explicitly written in appendix |X[ eq. (jHUj) . 



Results 

Figure shows r] as function of m\ and of m^/T^H- Although it has been obtained for 
mjvi = 10 10 GeV, other values of mjvi < 10 14 GeV would essentially give the same result. 

We see that the final baryon asymmetry is strongly suppressed if T RH m Nl . This 
is due to entropy release from inflaton decays, which gives a ~ (rn Nl /T R n) 5 suppression 
of rj. Furthermore during reheating the universe expands faster than during the standard 
thermal phase: H/ H staxidwA ps (T/T RH ) 2 - This makes both N x production and washout less 
efficient, increasing the value of fhi at which leptogenesis is maximally efficient, as can be 
seen in fig. [T^l 

Making use of the bound in eq. (j41)l [3T1| ITO! ITT] and marginalizing over m\ we can 
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Figure 13: Isocurves of the efficiency parameter r\ = lO" 6 ' -5 '---'^ 1 for leptogenesis with 
reheating as function of (mi, m^/Tnu) in the SM (left) and MSSM (right) for = 
10 10 GeV (the plot would be only slightly different for any m^ <C 10 14 GeV). 



therefore derive the bound on the reheating temperature shown in fig. This bound 
holds assuming that thermal leptogenesis generates the observed baryon asymmetry, that 
the inflaton decays into SM particles rather than into right-handed neutrinos, and that 
right-handed neutrinos are hierarchical, m^ <C TriN 23 - In the case of the SM, the lowest 
reheating temperature allowed for successful leptogenesis turns out to be 2.8 x 10 9 GeV, 
while in the case of the MSSM the lowest value is 1.6 x 10 9 GeV. However, there are 
various reasons to suspect that the reheating temperature is small in locally supersymmetric 
theories. Indeed, gravitinos (and other dangerous relics like moduli fields) are produced 
during reheating. Unless reheating is delayed, gravitinos will be overproduced, leading to 
a large undesired entropy production when they decay after big-bang nucleosynthesis [IE]. 
The limit from gravitino overproduction is Trh ^ 10 9 to 10 10 GeV, or even stronger |^9*] . 
This upper bound is at odds with the lower bound we have computed to achieve successful 
leptogenesis. Fig. Eh shows that this conflict can be avoided if N± and/or N± decayed while 
giving a substantial contribution to the total energy density of the universe and fhi ~ fh\. 
In such a situation leptogenesis depends on the precise value of the initial Ni, Ni abundancy, 
unless it is dominant. It can be realized if the inflaton decays dominantly into right-handed 
(s)neutrinos, or if Ni acquires a large vacuum expectation value, as discussed in the next 
subsection. An_ alternative solution to solve the gravitino problem, maintaining a thermal 
abundance of N±, is to rely on "soft leptogenesis" [T^ . 

We can further elaborate on the results presented in figures El El by making simple 
analytical approximations. Since we are interested in the effects of reheating, we consider 
the case m Nl > T RH , and since we are studying lower bounds on T RH , we restrict ourselves 
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to the most favorable case in which Tmax > n^Ni- 

The efficiency factor 7] receives three kinds of contributions, 

V = r/ab^eq^RH- (54) 

Here rj^ measures the Ni abundance before decay, relative to the equilibrium density. In 
order to estimate it, we first define 

K{T) = I (55) 

where Y is the N-y decay rate (r = (G Ffhimj^J / (2\^2n) at T < m Nl and 
r = (G j pmim^r i )/(2v / 2vrT) at T > m^) and H is the Hubble constant, with H = 
[57r 3 ^(T)]/[9^(T RH )] 1/2 T 4 /(T ] | H Mpi). Under the assumption 10 that the right-handed neu- 
trino density is much smaller than the equilibrium density , and taking the inverse 
decay as the dominant production process, in the relativistic limit we obtain jSH] 



d(n Nl /T 8 ) 8n 



cq 



dT 3 T 9 

Using K oc T~ 5 and n c ^ oc T 3 , we find 

n Nl 4 



K^r- (56) 



K. (57) 



Next, we define 

K. = K(m.,r) = 

3 x 10~ 3 eV \m Nl 

If ^> 1, the right-handed neutrinos reach the equilibrium density before they become 
non-relativistic and r\, A ^ = 1. If <C 1, from eq. (jBTj) we obtain 7? ab = (4/15)if*. 

The next coefficient in eq. (J5%jl is 77 eq , which measures the out-of-equilibrium condition at 
decay. If K* <C 1, the right-handed neutrino is decoupled when it becomes non-relativistic, 
and rjeq — l. If K* 1, ?7 C q can be estimated by computing the JVi density at the temper- 
ature Tf at which the processes that damp the baryon asymmetry go out of equilibrium, 

2£ (m^Y 2 IT , 



* = imvt) e ■ (59) 

We assume that the dominant process erasing the asymmetry is the inverse decay, with 
Tid = (m^/T) 3 / 2 exp(— mAr 1 /T)Gi?mim^ l /(8y / 7r). Then Tf is given by the condition 

T ID = H\ T=Tf . (60) 

If Tf > T RH , eq. (J5Uj) corresponds to A^(m A r 1 /Tj) /3 exp(—m Nl /Tf) ~ 1, with j3 = 11/2 
which, in the range 1 C < 10 4 is approximately solved by m Nl /Tf ~ a(lnfC) b , with 
a = 10 and 6 = 0.5. If Tj < T RH , the usual radiation-dominated epoch determines the 



30 



10 1 



> 10 11 
a 

s 

^ 10 10 



10 y 



10 y 



SM 



allowed 



10' 



10' 



niff, inGeV 



10' 



MSSM 




wijv inGeV 



Figure 14: Lower bounds on mjVi and Trh /rom leptogenesis in the SM and MSSM. 

dynamics and we find an analogous solution with (3 = 7/2, a ~ 5, b « 0.5, and -?C must be 
computed using the usual Hubble parameter H ccT 2 . 

Finally, ?]rh measures the dilution caused by the expansion during the reheating phase. 
Therefore rj RH = (T RH /m Nl ) 5 for < 1, 77 RH = (Trh/T/) 5 for A"„ > 1 and 7) > T RH , 
and ?7rh = 1 for if t > 1 and Tf < T RH . 

Putting together the different contributions to r], we obtain, for Tf > Trh 



77 « 20 



77 » 0.1 
10" 3 eV 
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Leptogenesis from infiaton sneutrino decays 

In supersymmetric seesaw models there is a distinctive possibility that infiaton itself is a 
scalar superpartner of the lightest heavy neutrino [oTJl 1ST] . This is an interesting scenario 
because both the reheating of the Universe and the thermal leptogenesis efficiency depend 
on a single neutrino parameter fh\. Therefore this is a predictive example of a realistic 
scenario of the early universe. In this case there is an additional source of the lepton 
asymmetry from the direct decays of the infiaton sneutrino. 

We do not study in detail how the sneutrino condensate decays, and assume a decay 
width r<£ = (T = 0) with CP asymmetry e\ = e^ {T = 0). This is not correct if 
2rh ?Z wiNi P^j - that, in our case, happens for mi>10 _3 eV. However if mi>10 _2 eV 

10 The expression we are using for T is not correct at high temperatures, where the Higgs decay is the 
relevant process. However, for this qualitative discussion, the approximation is adequate, since the right- 
handed production is dominated at temperatures T ~ m^, where we can take T = (GFrhim 3 Ni )/ (2^/2ttT). 
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Figure 15: Lower bounds on from sneutrino inflaton leptogenesis (left) and the effi- 
ciency of thermal leptogenesis in this scenario (right). 

thermalization is so efficient that details of the past thermal history negligibly affect our 
final result. 

In order to study this scenario we solved the Boltzmann equations for the N± = Nx ± 
N{ and N\ abundancies, and for the N\ condensate. In our calculation we take into 
account the temperature dependent interaction rates and CP asymmetries, and the on- 
shell resonances are correctly subtracted. The Hubble constant H and the parameter Z 
are the obvious supersymmetric extensions of eq.s ()80H82[) including the reheating effects 
of thermal (s)neutrinos. 

Assuming as before hierarchical light neutrinos and using the maximal CP asymmetry 
of eq. (jH| for supersymmetric case, the solutions to the Boltzmann equations are presented 
in fig. [T5l where we plot the lower bound on as a function of rh\ from the observed 
baryon asymmetry of the Universe. This parameter space has two distinct regions. The 
one denoted by A is the region of purely thermal leptogenesis. The corresponding curve in 
[S2] was obtained with wrong subtraction of on-shell resonances and with a constant CP 
asymmetry. For the region A we plot also the isocurves of the leptogenesis efficiency 77 
(right plot) which decreases very fast for small rh\. This is because of the suppression of 
Trh in that region of the parameter space. 

The region denoted by B is the one of direct leptogenesis from the inflaton decay |BT] . 

Between those two regions leptogenesis is a mixture of the two scenarios. Thus the 
inflaton sneutrino scenario allows to lower the lower bound on and on Trh from 

successful leptogenesis over a large region of rh\. This is very desirable from the point of 
view of the gravitino problem, as discussed in the previous subsection. 
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7 Conclusions 



We have performed a thorough study of thermal leptogenesis which, at present, is one 
of the most attractive mechanism to account for the baryon asymmetry of the universe. 
The final prediction of leptogenesis for the baryon asymmetry can be written in terms 
of the CP-asymmetry at zero temperature, e^, and of the efficiency t] of leptogenesis as 
n B /s = -1.37 10~ 3 e Nl r] in the SM, and as in eq. flUD for the MSSM. Figures El (SM) and[TT] 
(MSSM) show rj as function of the relevant unknown high energy parameters, i.e. the mass 
mjvi of the lightest right-handed neutrino Nx and fax, its contribution to light neutrino 
masses. All the new effects discussed in this paper have been added: cumulatively the final 
baryon asymmetry gets typically corrected by a order unity factor with respect to previous 
studies. For example, rig gets roughly doubled if fhx ~ (Am^^) 1 ' 2 . Since individual 
terms give larger corrections to the final result, in general, it is necessary to include these 
corrections to obtain a trustworthy result. The most important sources of corrections are 
summarized in section l4~Tl 

Although thermal leptogenesis allows to compute the baryon asymmetry in terms of 
particle physics, a few relevant parameters are presently unknown. Improving on this issue 
is as crucial as hard. In the meantime, by making some assumptions on the high-energy 
parameters (the most relevant one being that right-handed neutrinos are hierarchical) one 
can get interesting constraints [201221 • Including all the new effects discussed in this paper 
and combining uncertainties at 3a, we have found that successful leptogenesis needs 

( 2.4 x 10 9 GeV if Nx has zero 
mjvj > \ 4.9 x 10 8 GeV if Nx has thermal initial abundancy at T ^> 
( 1.7x 10 7 GeV if Nx has dominant 

and rris < 0.15 eV, where 7773 is the heaviest left-handed neutrino mass. In the MSSM we 
get similar values. 

Furthermore, in inflationary cosmologies, we obtained a lower bound on the reheating 
temperature, Trh > 2.8 x 10 9 GeV assuming that inflaton reheats SM particles but not 
directly right-handed neutrinos. Within the MSSM the bound is Trh > 1-6 x 10 9 GeV, 
which is at odds with the lower bound from gravitino overproduction. This seems to suggest 
that one has to rely on alternative (non-thermal) mechanisms to generate right-handed 
(s)neutrinos after inflation (like the sneutrino condensate studied at pageETJ), or to invoke 
leptogenesis with degenerate right-handed neutrinos 123122] or "soft leptogenesis" [12] (that 
we study at page l23|l . 

We stress that all these constraints are based on untested assumptions and therefore 
cannot be considered as absolute bounds. 

* * * 

Finally, we would like to emphasize some weak points and possible refinements of our 
analysis. At fh\ ^> 10~ 3 eV the relevant abundances are close to thermal equilibrium, 
suppressing the dependence on initial conditions. In this region we are not aware of any 
missing effect larger than ~ 10%. Our inclusion of thermal effects focussed on thermal 
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corrections to kinematics: by resumming corrections to propagators we could study effects 
which become large at T > m^. We approximately included corrections to couplings by 
renormalizing them at ~ 2nT. Although this is a significant improvement with respect to 
previous computations which used couplings renormalized at the weak scale, a somewhat 
different approach could give a more precise result valid for fh\ ^> 10~ 3 eV. As explained 
in the text, one should concentrate on computing corrections to the N\ — > HL decay rate 
at relatively small temperature, T < m^, without making our simplifying approximations, 
with the goal of including all few % corrections of relative order ~ g 2 / it 2 and \ 2 /ir 2 : AL = 1 
scatterings and their CP-asymmetry, three-body decays Ni — > LQ3U3, LHA and radiative 
corrections to N\ — > LH decay and its CP-asymmetry. 

At fh\ ;$ 10~ 3 eV the final result depends on initial conditions. Starting with zero initial 
abundancy, the final baryon asymemtry also depends strongly on thermal corrections to 
the CP-asymmetry e^. Unfortunately we found that, as the temperature rises, thermal 
corrections first reduce, then enhance, reduce and finally enhance e^. Since the correction 
does not go in a clear direction, a more accurate computation might be welcome. 
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Hambye, M. Laine, M. Mangano, M. Moretti, M. Papucci, M. Passera, M. Pliimacher, R. 
Rattazzi, F. Vissani. 



A Boltzmann equations 

Elastic scatterings keep the SM particles in kinetic equilibrium, so that their energy dis- 
tribution approximatively follows the Maxwell-Boltzmann distribution / = e~ E l T (unless 
otherwise specified we neglect the slightly different energy distributions of bosons and 
fermions. This is a good approximation because / < 0.05 at the average energy (E) ~ 3T). 
Each particle species is simply characterized by its total abundancy, that can be varied only 
by inelastic processes. When they are sufficiently fast to maintain also chemical equilibrium, 
the total number n cq and energy density p cq of ultra-relativistic particles at temperature T 
are 



(2tt) 3 j tt 2 ' r a J (2tt) 3 j it 2 

where g is the number of spin degrees of freedom. All SM particles have g$M = 118, and a 
right handed neutrino has = 2. 

In kinetic equilibrium, the phase space density is fn/n cq . The Boltzmann equation 
describing the evolution of the total abundancy of a particle N is 

h N + 3Hn N = -^2[Na---++ij---], (62) 



where H = a/a = y^8-Rp/3Mj 



pd 



n^n a ■ ■ ■ Tij rij ■ ■ ■ 
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A symmetry factor should be added when there are identical particles in the final state or 
in the initial state. 7 eq is the space-time density of scatterings in thermal equilibrium of 
the various interactions in which a N particle takes part: 

7 eq (A r a ^ ij) = J d £ Nd £ a f N f a J d p. d $. (2ir) 4 8 4 (P N + P a -P t - P,)\A\ 2 , (63) 

where dpx = d 3 px/2Ex(2n) 3 and \A\ 2 is the squared transition amplitude summed over 
initial and final spins. Neglecting CP-violating effects, the inverse processes have the same 
reaction densities. 

Having neglected Pauli-blocking and stimulated emission factors, and assuming that 
| A | 2 does not depend on the relative motion of particles with respect to the plasma, the 
expression for the scattering rates 7 eq can be conveniently simplified. 

For a decay, eq. (J55|) reduces to 

7 eq (iV - ij ■ ■ ■ ) = 7 cq (u ••■-#)= ^|44 Far, (64) 

where z = m^/T and IV is the decay width in the rest system of the decaying particle. 
The Bessel K functions arise from the thermal average of the time dilatation factor m^/E. 

The thermal rate 7 cq of a two-body scattering can be conveniently rewritten by multi- 
plying eq. fU) times 1 = / d 4 P J 5\P - P N - P a ) 

^ {Na «_ ij . . . ) = j_ r £l &{Na _ i:j . . e -Po/T, (65) 

87T J (27TJ 4 

Here a = 2s A[l, m%/s, m 2 /s] a is the 'reduced cross section', a is the total cross section 
summed over final and initial quantum numbers (spin, gauge multiplicity, . . . ), X[a, b, c] = 
(a — b — c) 2 — Abe. If the total cross section a depends only on s (and not on the thermal 
motion with respect to the plasma) we can use d A P = lix^/Pl - s ds dP to get 

7 eq (iVa - ij ■ ■ ■ ) = r ds s 1 ' 2 a{s) K x (^) , (66) 

where s is the squared center of mass energy, s min = max[(m A r + m a ) 2 , (m, + rrij + ■ ■ ■ ) 2 ]. 

The 3H term in the Boltzmann equations, eq. (|62j). accounts for the dilution due to 
the overall expansion of the universe. It is convenient to reabsorb it by normalizing the 
number density n to the entropy density s. Therefore we study the evolution of Y = n/s 
as function of z = m^/T in place of time t (H dt = dlnR = dlnz since during adiabatic 
expansion sR 3 is constant, i.e. R oc 1/T). The Boltzmann equations become 

zHs^ = - l Na ~ ( 67 ) 

a,i,j,... 
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Leptogenesis 



We now specialize to leptogenesis. Neglecting sphalerons, the scattering processes relevant 
for leptogenesis are: 11 



where 



AL 
AL 



A. 



1 : D =[Nx*r-> LH] S s = H s + A s S t 

2 : N S = [LH <-> LH) N t = [LL <-> M] 



[LiVi «-> 



Q 3 U 3 ], 2H t = [N 1 U 3 ^Q 3 L} + [N 1 Q 3 ^U 3 L}, 
HA], 2A t = [N 1 H ^ AL] + [N 1 A^ HL}. 



(68) 



(69) 



We separated AL = 1 scatterings S Sjt into top (H s>t ) and gauge contributions (As )t ). The 
relative Feynman diagrams are plotted in fig. H] and computed, including finite temperature 
effects, in appendices iBl and IU1 

We assume that Ni, L and L can be out of thermal equilibrium, while Yx 
X = {H, Q 3 , U 3 }. The Boltzmann equations are 



Y x q for 



zHsY^ 
zHsYi 
zHsY[ 



D 
D 



-D-D 

N„ - 



Ns 



S s S s St 
N t -S s + S~ t 
W t -S~ s + S t . 



We write the decay rates in terms of the CP-conserving total decay rate 7d and of the 
CP- asymmetry <C 1: 



so that 



D 



2 



Y, 



Ni 



7 oq (iVi 



LH) = Y«(LH -> NO = (1 + ejVl ) 7 D/2, 
L#) = 7 cq (^ -> iVi) = (1 - e Nl )j D /2 



(70) 



A, 



7d 
2 



Here Y^*, Y^ q and Y^ q are fermionic equilibrium densities each with 2 degrees of freedom, 
and therefore are all equal to high temperature. Keeping only decays and inverse decays, a 
baryon asymmetry would be generated even in thermal equilibrium, since CPT invariance 
implies that if Ni decays preferentially produce L, than inverse decays preferentially destroy 
L (or, in formulae, D — D ^ 0) [SI]. In order to obtain Boltzmann equations with the correct 
behavior one needs to correctly address some subtlety, as discussed in ref. jSH] in the context 
of baryogenesis. 

11 We add AL = 1 scatterings involving gauge bosons. We neglect three-body decays N\ — ► LQ3U3, LHA 
and radiative corrections to N\ — > LH decay and its CP-asymmetry, although in most of the parameter 
space they enter at the same order as AL = 1 scatterings, giving ~ g 2 /^ 2 ~ few % corrections. Ref. [32] 
suggests that AL = 1 scatterings have, at T = 0, the same CP-asymmetry as decays. If true this gives % 
corrections, that we prefer not to include. 
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Subtractions of on-shell propagators 



At leading order in the couplings, 2^2 scatterings must be computed at tree level and 
are consequently CP-conserving. However N s , the LH <-» LH scattering rate mediated 
by s-channel exchange of N\ shown in fig. EJl, must be computed by subtracting the CP- 
violating contribution due to on-shell Ni exchange, because in the Boltzmann equations 
this effect is already taken into account by successive decays, LH <-» N <->• LH. Since 
the on-shell contribution is ^~ shcl \LH -> LH) = ^(LH -> iVi)BR(iVi -> Z#), where 
5i2(iVi -> L#) = (1 - ejvJ/2, we obtain 

^(LH^LH) = lNs -(l- £Nl ) 2 lD /4, (71) 
^{LH-tLH) = lNs -(l + e Nl ) 2 lD /4, (72) 



so that 



iV s = ^^(LH ^ LH) - ^^(LH ^ LH) (73) 
^-^U^d + ^U (74) 



— ySq \^ Ws ^- y T cjvi CD T ^V c 7Vi, 

L 

having defined the lepton number asymmetry Yc = Yi — Yi and used Yj, + Yi — 2Y^ q + 
O(ejvi). Expanding at leading order in gives the Boltzmann equations 



Y, 



zHsY^ = -(^-l)( 7D + 2 75s + 4 75t ) (75) 



zHsY' c = 7D 



Y 



yv, a y, 



V eq / OV eq 



^(2 7 - b + 2 75t + 7Ss ^), (76) 



where 7 ™ b = 7A r t + 7A r s - 7d /4. 

Alternatively, one can simply not include the decay contribution to washout of Yc 
because it is already accounted by resonant decays. Then one gets 

zHsY'c = lD e Nl - lj - ^ ^2 lN + 2 75t + Issyj^j , (77) 

which is equivalent to eq. (|7T)j) . Therefore it is not necessary to compute subtracted rates. 
For our purposes it is convenient to compute 7 at as 7 ™ b + 7r>/4. Since 7j d is a simple and 
important quantity, we will compute it accurately going beyond the Boltzmann approxi- 
mation. We now discuss how to directly compute 7™ b , and show that it is subdominant, 
unless the Ni Yukawa couplings are large. 

We can directly compute 7 ™ b = 7 tvs — 7z?/4 by subtracting to the intermediate N\ 
propagator the resonant part in the narrow-width approximation 12 

\D Nl \ 2 - \Dt\ 2 = \D Nl \ 2 ^— 6{a - m 2 Ni ) D N = g-L ( 78 ) 

m Nl T Nl s - m% + im Nl T Nl 



12 In numerical computations one employs any representation of S(x) that, like 5(x) = (2e 3 / V) / \x 2 + e 2 ) 2 
and unlike 5(x) = (e/ir) / (x 2 + e 2 ), goes to zero faster than the propagator away from the peak. The value of 
e can be conveniently chosen to be somewhat smaller than the width of N\ , although this is not necessary 
if it is narrow. In this limit, which covers almost all the parameter space where thermal leptogenesis 
can generate the observed baryon asymmetry, one can simply set T^/m^ = e, getting the subtracted 
propagator of eq. 1)88(1 . and assign to the width any sufficiently small value. 



37 



as discussed, in a different context in ref. [SH]. Using eqs. (J70*Jl and and recalling that, 
for s-channel exchange a(LH — > Z/f) = 87rr^ \D^ l | 2 , we obtain 



f°d a s^a(s) Kl (^) 

Smin 



sub _ T f°° j„ „l/2 T ^ 



rjViJTljvJ-DiVi 



5(s -m^J 



7at s -^. (79) 



This corresponds to the result previously obtained. 
The subtraction employed in ref.s EH E2I is 



\D Nl \ 2 - iRe^J 2 = I^J 2 - (ImD^) 2 ^J 2 - 1^— 5( s - m ^ ) 

2 m Nl L Nl 

which leads to Yns ~ 7tvs — 7r>/8. This corresponds to subtracting only 1/2 of the on-shell 
contribution, thereby over-estimating washout by 3/2 when 7^ is the dominant process (i.e. 
when the neutrino Yukawa coupling is small). As shown in fig. the properly subtracted 
rate has no spurious peaks around the resonance region, in contrast with the result of 

refs. gaum. 



Sphalerons 

Finally, we have to include sphaleronic scatterings, that convert the lepton asymmetry 
into a baryon asymmetry. This is conveniently done by converting the Boltzmann equa- 
tion for Yjc into a Boltzmann equation for Yb-c'- since Yb-c is 110 1 affected by sphalerons 
we only need to find the relation between Yb-c and Yc- At temperatures larger than 
10 10 GeV sphaleronic scatterings are expected to be negligibly slow with respect to the 
expansion rate of the universe, so that Yc = —Yb-c- At lower temperatures sphaleronic 
processes keep thermal equilibrium and the relation would become Yc = — (63/79) Yb-c 
(Yc = — (9/8)Y#_c) if all SM Yukawa couplings were large (negligible). In reality some 
couplings mediate equilibrium reactions (y t , . . .) and some others are negligible (y e , . . .) so 
that without making approximations one cannot ignore flavour and must proceed as in 
ref. jH]. In particular we stress that in order to study how the three generations share the 
lepton asymmetry one must consider the evolution of the full flavour 3x3 density matrix. 
Within 10% accuracy we may approximate Yc ~ — Yb-c an d solve the Boltzmann equation 
of eq. (|80c|) . After computing Y B -c the baryon asymmetry Y B is obtained by means of 
eq. flUID. 



Inflaton and N\ reheating 

We now add one refinement. We described in section El how the Boltzmann equations are 
modified by the presence of a field 0, whose decays into SM particles reheat the universe. 
Proceeding along the same lines we also take into account that reheating due to N\ decays 
might be non-negligible. Terms of relative order Pn 1 /pr are neglected in the 'standard' 
Boltzmann equations: in thermal equilibrium this factor is small, P^/pr^ 9Ni/ 9* ~ 0.02 
since Ni is one out of many more SM particles. Away from thermal equilibrium it can be 
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sizable. Including these effects, the Boltzmann equations become 



HZz—^- = SHpj, - r^p^ , (80a) 
sHZz^- = 3sH(Z-l)Y Nl -(^-iy iD + 2 1Ss + A 7st ) (80b) 



sHZz—^ = 3sH(Z-l)Y B _ c - lD e Nl i^-l 



■JVi 

Y B~C (lD sub „ ^JVi i 



where s and p# are the entropy and energy density of SM particles, 



is the Hubble constant at temperature T and 



Id + ^Iss + ^Ist Pn ± ( Y i 



4Hp R AH p R n e l \Y£ 



(82) 



takes into account reheating. Since interactions with SM field can keep iVi close 
to thermal equilibrium, its contribution to Z is not just given by IT 'n 1 Pn 1 / '4i? Pr = 
(ln/4Hp SM ) {p* /„£ ) (Y N JY£) . 

In the next appendices we will see how thermal effects can be included by modifying 
the scattering rates and the CP-asymmetries, but not the form of the equation themselves. 



B Thermal corrections to decays and scatterings 

We present the temperature-dependent decay rates and cross sections that generate and 
washout the lepton asymmetry. Since we consider temperatures T ~ much above the 
electroweak scale, the thermal masses of the Higgs doublet, lepton doublet, top quarks and 



electroweak gauge bosons are j 


pven by [2] 








m H 


3 9 1 2 1 2 1 * 


(83a) 


m| 


- 3 2 ^ 2 + 3 2 ^, 




(83b) 


22 


1 2 

= ~ 6 9i + 


3 2, 

32^ 2 + 


1 2 1 2 

288^ y+ 16 yt ' 


(83c) 


m u 3 


1 2 

= 


1 2 

18^ + 


1 2 


(83d) 


m W 
J>2 


11 2 
- 12 ^2» 


m| 

J>2 


11 2 

- 12 ^y, 


(83e) 
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where all couplings are renormalized at the RGE scale \x = 2nT. The Higgs boson, lepton 
and quark masses are functions of T, even if not explicitly denoted. We neglected their 
zero-temperature values. At leading order the quartic Higgs coupling A is given in terms of 
the zero-temperature Higgs mass nih as X(fi — v) = (rrih/2v) 2 where v = 174 GeV. The top 
Yukawa couplings is similarly given by yt(fi = v) = m t /v. Their RGE running is known up 
to next-to- leading order [H7j . 

In the MSSM the relevant thermal masses are |2Sj 



■in 



7^2 

2 
L 



111 



J^2 
J^2 



2 

m H 
O rtu 

2 

m T 
2— 

2 

Ji2 



m 



u 3 



J^2 



^ 2 ^ 2 

8<? 2 + - 8 9y ~ 

3 o 1 O 

8 92 + 8^ y ' 

2 o 3 9 

3^3 + g^2 4 



3 A 2 

4 A*, 



1 2 ! X 2 

72^ y + 4 A *' 



(84) 



3^ + ~/ Y + 2 A? 



Here and in the following appendix the neutrino couplings Y v are renormalized at the high- 
scale. The one-loop RGE equations for the Majorana neutrino mass matrix valid from 
the Fermi scale (below which it is not affected by quantum corrections) up to are 



(4tt) 



j dm 
din// 



m(A — 3#2 + 6A 



m.{YlY E ) T + {YlY E ). 



in 



55) 



in the SM and 



(4tt) 



, dm 
dln/x 



m 



2^ _ 6^ 2 2 + 6A 2 ) + m • (yj • Y E ) T + (Y] • Y E ) 



m 



56) 



in the MSSM. Higher order effects (two-loop RGE, thresholds) are partially included in our 
codes. Here Y E and A are the charged lepton and Higgs couplings. The solution of these 
RGEs is described in section 13.11 The Yukawa couplings of right-handed neutrinos N\ give 
extra RGE effects at scales above m^. We neglect these effects, as large Yukawa couplings 
anyhow lead to an exponentially small efficiency for leptogenesis. 

The Ni mass does not receive thermal corrections, as long as we neglect the relevant 
neutrino Yukawa couplings, which are indeed small in most of the interesting parameter 
region, since (Y^l 2 < 3 x KT^m^/lO 10 GeV)(mi/l(T 3 eV). We define 



m Nl 
T 



x 



m 



A'l 



&H,L,Q,U,W,B 



m H,L,Q 3 ,U 3 ,W,B{ z ) 



111 



F 2 



m 



57) 



D 



A'l 



X 



ia 



1/2 ' 



and A[a, b, c] = (a — b — c) 2 — 46c. 



I /-)2 I sub 
l-^JVil 



[X 



[(x 
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Decays 



At low temperature > m H + m L so that one has the usual Ni decay with total width 
T Nl = ^- (Y^Y^) 11 A^[l, a H , a L ](l - a H + a L ) m Nl . (89) 

S7T 



Analytical approximate solutions of the Boltzmann equations |2II| show that in the most 
interesting region (fhi ^> eV and small Ni Yukawa couplings), the efficiency of thermal 
leptogenesis is mainly controlled by 7^ = 7™ b + 7d/4 — 7d/4. Therefore we compute the 
thermally averaged N\ — > LH decay rate 7^ accurately. Instead of using the Boltzmann 
approximation of eq. ()64|). we compute 7^ using the full Bose-Einstein fu and Fermi-Dirac 
fi, /jy distributions, including stimulated emission and Pauli blocking factors. The relation 
7jv — 7d/4 remains unaltered where now 

Id = 2jdp N dp L dp H f L fHfNe EN/T (2nY5\P N -P L -P H )\A\ 2 (90) 
= dE N m Nl JE 2 N -m 2 N J N e E ^ T - dcosB f L {E L )f H {E H ) 

n Jm Nl V Z J-l 

where dpx = d 3 px/2E x (2,7r) 3 is the relativistic phase space, 9 , E^,p L are the decay angle, 
the L energy, the L momentum in the N± rest frame and El = ^E\ + \/"f 2 — lp° L cos 6*0 
(7 = E N /m Nl ). The integral in dcos9 can be done analytically. 

As expected, the Boltzmann approximation is accurate at low T <§.lrriN 1 and differs 
from the full result by a few 10% at T ~ tyin^ ■ Actually, for the specific SM values of the 
thermal masses, the Boltzmann approximation is accurate within 10%. 

We stress that our inclusion of thermal effects is only approximate. In the present work 
we focussed on those corrections which become large at T> m^. 

In particular, at sufficiently high temperature the Ni decay becomes kinematically for- 
bidden at higher temperature by the H thermal mass. When m# > + tul, the 
H — > LNi decay is allowed. Its width is 



(y v yJ)nAa[aj,l,a £ ]— 5 — ~ m H- (91) 



H 



Despite the change in the decay process the Boltzmann equations keep the same form, with 
the N\ decay rate replaced by the H decay rate. 

Even including thermal effects, at intermediate temperatures all 1 «-> 2 and 3^0 
processes are kinematically forbidden, so that 7^ = 0. 



AL = 2 scatterings 

We now consider lepton-number violating scatterings: the AL = 2 processes LH^ <-» LH 
and LL <-» H^H (middle row of fig. EJ). 

As explained in appendix [X] we compute the LH «-> LW subtracting the contribution 
from on-shell intermediate N\, to avoid double counting with two-body decays. At low 
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temperatures the computation is similar to the T = case: Ni can be on-shell in the 
s-channel diagram of fig. EJi so that one must replace the s-channel propagator = 
l/(s — m 2 N + im^T^) with a subtracted propagator D snh . 

At high temperatures, when H — > N\L replace N\ — > HL decays, the situation becomes 
more tricky. JVi can be on-shell in the u-channel diagram of fig. Et. As this corresponds 
to on-shell N\ in the decay H — > N±L, it has to be subtracted similarly to the s-channel 
resonance. The imaginary part of the N\ propagator which renders finite the LH <-» LH^ 
rate is no longer given by Ni decay, but by thermal absorption of Ni, given by (SHj 



-Imn Wl (E 



— -^jvTjvi (-S'iv) — - 



dp L #W (2tt) a 5\P n + Pl - Ph)\A\Vl + fn] 



Unlike a decay at T = 0, for which F^^En) = ^N 1 { m N 1 ) m N 1 / En, the Ni width at finite 
temperature depends on thermal motion of N% with respect to the plasma, giving rise to 
lengthy expressions. For simplicity we give the expression corresponding to the narrow- 
resonance limit IVl — > 0, which is always valid. In fact, if the N% Yukawa coupling is large, 
Ni gets a thermal mass avoiding H — > N\L decays. 

The reduced LH «-> LH^ cross section is given by 



~ sub 



V 111 



Atix 



1 + (1 - 2a H + x)(ReD Nl - 3f)) In 



R 2 /x 2 



+^(fl 2 H + (a L - x) 2 - 2a H {a L + x)) |D 



x- 



)2 I sub r>2 
JVi I K 2 



Rj + e 
+ 2xReD Nl 



+ 



(92) 



2(1 + x- 2a H )R 2 Ri 

[R 2 /x 2 + e][R 2 + e] 



3£(2x 



[x 



a-H + a L ) 2 (ReD Nl - f) 



where R\ = 1 — 2a h — 2a^ + x and i?2 
For the AL = 2 scattering LL 

(Y u Yj) 2 u x-2a L 



x - 



(an — o-l) 2 and e is any small number, e <C a\. 



0~Nt 



2tt 



x 



Ri 



HH we obtain 

r 3.9-1 

L _ r £ 2 _L_ ( 



30 In 



Pi 



7?i + 1 - r 



(93) 



where r = yj(x — Aa H )(x — Aa L ). The parameter £ takes into account scatterings mediated 
by heavier right-handed neutrinos m N2 3 ^ m Nl , as explained in section HJ 



AL = 1 scatterings 

Let us now consider the AL = 1 processes LN QaL^ and C^iV <-> Q3L (bottom row of 
fig. EJ). We can neglect the small difference between the thermal masses of Q3 and C/ 3 (see 
fig. QJi), setting ac/ ~ oq. The LiV (53P3 reduced cross section is 

»h. = ^-(Y^) n y 2 t { - - 1 7 ^ ~ 2QQ) ^[(1 + a L - ^ - 4a L ][l - 4a Q /x], (94) 
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The U%N *-> Q3L and the Q3N U3L cross sections are equal and given by 

1 1 



t + — i_ — (1 — an + a>L){ a H — 2a 



Q) 



'1 - 2a H + a L + 2a Q ) In 



t + — an 
t- — an 



(95) 



where 



1 

2x 



a Q + x - (a Q - x) 2 + a L (x + a Q - 1) + 



±\/[ a Q + 0* ~ !) 2 - 2a Q (l + x)][al + (x - a Q ) 2 - 2a L (x + a Q )} 



(96) 



Neglecting thermal masses, all these interaction rates agree with those used in the 
literature (up to typos present in older papers). As explained in the text, we must include 
AL = 1 scatterings involving gauge bosons. We do not compute the full thermally corrected 
rates, since we should take into account thermal motion of A and L with respect to the 
plasma, which gives a complicated result. We cannot fully neglect thermal effects, since 
exchange of massless H and L would give an IR divergent result. Therefore we keep 
thermal masses of A, L, H only when they act as regulators of IR enhanced contributions. 
This approximation is accurate at T m Nl . At larger T it neglects terms suppressed by 



higher powers of g 2 
a(NiL — > HA) = 



1/3. The result is 
3^(^)11 



167TX 2 



2t(x -2) + (2-2x + x 2 ) ln[(a L - t) 2 + e] + 
x(cll — t)(a,L + (Ilx — aw) + e(2 — 2x + x 2 



(97) 



+2- 



(a L -t) 2 + e 



a(LH -> NxA) 



a (LA -> NiH) 



3p 2 2 (^) 
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8irx(l — x) 

Zg 2 (YX)ii 
IQtxx 2 



2zln(* - a H ) - (1 



x 



ln(t 



t- 
-\t+ 



x 



CL\y — a H 1 



t 2 + 2t(x - 2) - 4(x - 1) ln(t - a H ) 



x- 



a w — 4a H 
an — t 



t- 

"1*4 



We wrote only the SU(2)x contribution. One must add the U(l)y contribution, obtained 
by substituting aw — > ay and 3 



1^2 l9 Y - 13 The expression [f(t)}\ + _ denotes f(t + /m 2 Ni ) - 



f(t_/m 2 Nl ) where t± are the usual kinematical ranges for t 
12 — > 34 scatterings: 



(^1 



in the various 



t± 



ml — m% — m 2 + m 2 ^ 2 



is 



m\ — m^) 2 



Is 



ml ± 



m\ — m 2 ) 2 



is 



13 Gauge scatterings have been estimated in a recent paper |82j . by introducing some infra-red cutoff, 
which should give a qualitatively correct result at low temperature. We can only compare the ratio between 
SU(2)i and U(l)y contributions, which is different from our value. We do not use simplified expressions 
for t±, valid when Tn\ WH <C s,m 2 Ni , because they not even accurate at low temperature, where small 



i Nl ~ m L W _ 



H is a relevant kinematical range. Due to the 1 — x factor, at T <C the LH — > NiA 



interaction rate is ^At ~ (ff/ 7r ) 2 7-D? of the same order as one loop corrections to the decay rate (that we do 
not include). 
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In some parameter range the process N\L — > HA can have on-shell L in the t-channel, that 
we have subtracted: in such a case one should use any finite value e <C a 2 L . Otherwise one 
can set e = 0. 

The rates that enter in Boltzmann equations are given by 

a As = a(LNi HA), a At = ~ [a{AL -> N,H) + a{HL iV^)] 



Resonances in s and it channels 

We here explicitly verify that the relation between decay and resonant scatterings, 
^on-sheii __ r^ D j\ ( U p to CP- violating corrections), remains valid without approximating 
thermal distributions with Boltzmann statistics, and that it applies to both s-channel as 
well as w-channel resonances. This last issue is non-trivial, as computing a cross section 
mediated by an unstable particle which can be on-shell is a difficult problem even at zero 
temperature and at tree level. In that case, the beam size and the width of the external 
unstable particle act gulator of the divergence [SU] • 

We first consider the s-channel case. The Ni <-> LH rates are 
7 eq (iV~i ^LH) = J dp Nl dp L dp H fmi 1 ~ h)(l + M (2tt) 4 5\P Ni - P L - P H )\A(N t LH)f 
7 cq (Lif - Ni) = J dp Nl dp L dp H f L f H (l - f Nl ) (27r) i 5\P Nl -P L - P H )\A{LH - N x )\ 2 . 

Using E Nl =E H + E L , we obtain f Nl (l-_f L )(l + f H ) = f L f H (l - f Nl ) = f L f Nl f H e E ^/ T . 
The on-shell contribution to the LH — > LH equilibrium interaction rate is 



' on— shell 



(LH - L#) = J dp H dp L dp L dp R f L f H (l - f L )(l + f s )\A(LH - m) 
ttS(s - m 2 Nl )' 



\A(N, - LH)\\27if5\P L + P H -P L - P R ) 



where the width that cutoffs the resonance is jSHl; the damping rate at finite temper- 
ature. We can rewrite the product of thermal distributions as: 

(1 - f L )(l + f s ) = f Nl e E ^ T [(1 - f L )(l + f s ) + f L fs\ (98) 

and insert 1 = J d 4 PN 1 S 4 (PN 1 — Pl — Ph), obtaining: 

f L fHf Nl e E ^ T \A{LH -> iV 1 )| 2 (2vr)V(P L + P H - P Nl ) ( 1 ) (99) 

[ dp L dp B [(1 - /z)(l + f B ) + hfs] \MNi - LH)\ 2 (2-k) a 5 4 (P Ni -P- l - P S ) 
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The integrals over final-state particles reconstruct the thermal width of N\ 

T% = — [ dp L dp R [(1 - f L )(l + f H ) + f L f H ] \A tiee \ 2 (2n)W(P Nl -P L - P B ) (100) 
fn Nl J 



In Boltzmann approximation the term in square brackets can be approximated with 1, and 
reduces to the standard expression, eq. (JHHJ). The integrals over initial-state particles 
reconstruct 7^, giving the desired relation 7™ s 611 = 7d/4. 

In an analogous way we can deal with u-channel resonance, present in LH — > LH 
scatterings at high temperatures when H decays to N\L while N\ is stable. At finite T 
the propagator of a particle involved in a 1 «-> 2 process gets an imaginary part, even if it 
is not the decaying particle. This thermal cutoffs the w-channel resonance and gives 
consistent Boltzmann equations (no £ asymmetry generated in thermal equilibrium). To 
show this fact we follow the same procedure. The interaction rates for H *-* N±L are 

7° q (# N ± L) = J dp n dp L dp Nl f s (l - f Nl )(l - f L ) (2ir) A 5\P s - P L - P Nl )\A(H LN{)\ 2 

7 eq (A r lL _^ S ) = J dp Nl dp L dp R f L f Nl (1 + f B ) (27r) 4 5 4 (P R -P l - P Nl )\A(N!L - H)\ 2 . 

Since E s = E Nl + E L one can show that fail - /j\r 1 )(l - h) = IlInA 1 + fa), so that, 
again 

^(H NlL) = y*( Ni l - i?) = ^(1 + e H ) (101) 
7 C( i(iV 1 L ^H) = 7 ^{H - N,L) = ^(1 - e H ) (102) 

where now 

lo = 2j dp n dp L dp Nl f n {l-f Nl ){l-f L ){2Ti)H\P fI -P L -P Nl )\A t UH ^ N^ 2 . (103) 
The u-channel on-shell contribution to LH — > LH is 

7 C n-shcii(^^ - LH) = f dp H dp L dp L dp R f L f H {l - f L )(l + f R ) 

\A(H - Nl L)\ 2 l %™ Nl \ \A(N lL - H)\ 2 (2.) H\P L + P H - P L - P s ). 
The product of the distributions can be rewritten, making use of Eg = E^ 1 + El as 

/l(i + f a ) = (1 - fm) [(/x)(l + fa) + fail - h)] (105) 

Inserting 1 = J d 4 P Nl 5 4 (P H — Pi — -P/vJ we get 

12-sULH - = J dpndpi {j0^^5{P 2 Ni - m 2 N 

f H (l - f Nl )(l - f L )\A(H - iV 1 L)| 2 (27r) 4 5 4 (P H - P L - P Nl ) ( nrn : ) 



2£ 7Vl r t 7V h l (£ 7Vl ) 

«^[(/l)(1 + /ff) + fail ~ MMNiL - H)\ 2 (2n)^(P Nl +Pl-Pr 
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Figure 16: Im[lQ(I{ + in N\ — > LH decay. The momenta flow in the direction in 
which the labels are written. 

The last terms reconstructs the thermal damping rate of JVi jSHj (when the decaying particle 
is not N\, but H), due to interactions of N\ with the L present in the plasma producing a 
H (and the inverse process), 



^- • 2 J dp L dp R [(f L )(l + f B ) + f B (l - f L )\ 
\A trcc (H -> LN 1 )\ 2 (2tt)^(P Ni + P l - P b ). 



(107) 



C Thermal correction to SM CP-asymmetries 

In this appendix we show explicitly the calculation of the CP asymmetry parameter e for 
both the Ni decay and the H decay, in the SM. Since the final result for leptogenesis turns 
out to be very weakly dependent on e#, we compute e# neglecting the motion of H with 
respect to the plasma. On the contrary we make the analogous approximation for ejvi only 
in the analytic expression presented in the main text, eq. (|30jl . and present here the full 
result, averaged over thermal N\ motion. 

We choose to work in the rest system of the plasma. The CP-asymmetry in i —> / 
decays which enters into the Boltzmann equations is (see eq. (|7D]l ) 



7 cq (^/)-7 cq (^7) _ Jd* Pi f(E t ) fd<j>e A \A\ 



7 eq(, _ /} + 7 e q( ^ _ /) J ^. J # \ A \2 



(108) 



where d 3 pi f{Ei) is the thermal distribution of the initial state, 
phase-space for i — > f decay, 

|A(i^/)| 2 -L4(W)| 2 



\A(i^f)\* + \A(-i^f)[' 



• is the standard relativistic 



(109) 



and A is the amplitude for the process computed at given quadri- momenta Pij. We define 
Ao{i —>■ f) as the amplitude of the tree level process and parametrize it as A = Xolo, 
where Ao contains the coupling constants. In the same way we define the amplitude up to 
the one-loop level as A = X I + With this notations e A at one-loop is expressed as 



Ao-^0 " 




_ Aq/ - 


hXth | 2 


Ao-^0 " 




+ |A57o - 


-Kh\ 2 



lAnl 2 



Unl 2 



(110) 
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In our case ii is the sum of two diagrams: a "vertex" correction and a "wave" correction 
(fig.sll6 | H7J ) . For both diagrams the couplings and the \A (i — > f)\ 2 factor are the same: 

o Im[(yty%] lm^h] 
e - 4 = " 2 (YV) n 2P N .P L 

where the quadri- momenta P/v and Pl are defined in fig. EH for Ni decay and in fig. IT7I for 
H decay. We have to compute Im[Io/i]. Using the cutting rules, we obtain 



lm[r o h] = ij J i- ( 112 ) 

cuttings 

Ii is the sum of two diagrams: the "vertex" and the "wave" one: I\ = l\ + 1^ . 



Computation of 

We consider here the case of Ni — > LH. In principle there are three possible cuttings 
(or six circlings in the notation of ref. [H2 121]) for l\ , and one (or two circlings) for 1^ , 
and none of them is forbidden at finite T. In fact, at finite T, energy is no more forced 
to flow from uncircled to circled vertices (see ref. [2Z|)- The reason is that, while the cut 
propagators at zero T are proportional to 0(E), at finite T they have a new contribution, 
see eqs. (fTU^ — f|^3|) . proportional to fp (or f B ) if the cut particle is a fermion (or a boson). 
This accounts for particles in the thermal bath which do not have the 0(E) function, since 
they can be emitted by the bath (positive energy) or absorbed from the bath (negative 
energy). However, recalling that we are working under the assumption that the Nj (with 
j ^ 1) are very heavy, the cuttings which involve the Nj can be neglected since they are 
exponentially suppressed by ttin /T. Moreover the graph with the iVi circulating in the 
loop does not contribute to the CP asymmetry since its Yukawa couplings are real (see 
eq. (jlll)) with j = 1). So, the only relevant cutting in l\ is the one with lepton line and 
Higgs line cut (fig. HBJ), as in the zero-temperature case. 

We compute here the cutting in the vertex part 

Im[I* lY] = j % J 0^D N (P L -P N +K) [D + H (P N - K)Dt(K) + D~ H (P N - K)D~ L (K)] T(K), 

(113) 

which is the first contribution shown in fig. (the single cutting in the figure stands 
for two possible circlings of the vertices). Here T(K) is the result of the traces over the 
spinor indices, D are the propagators (without numerators): is the propagator of the 
Nj (which we choose as the zero-temperature one, since we neglect Nj interactions with 
the plasma), and are the finite-T cut propagators of the Higgs and of the lepton 
respectively, see eq. (fTUj) and (J251) . 

1 

(P L -P N + Kf - m 2 Nj 

2tt5((P n - Kf - m 2 H ) [0(±(E N - u)) + f B (\E N - u\)} ,(114) 
2Tr6(6 L (u,k)) [0(±u) - f F (\u\)] , 



D N (P L -P N + K) = 

D±(P N -K) = 
Df(K) = 
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where 

5 L = [(l + a)uj + b] 2 -(l+a) 2 k 2 (115) 

and a, b are defined in eq. ([16|h Keeping only relative angles relevant for our computation, 
we can conveniently parameterize the quadri-vectors as 

U= (1,0,0,0), P N = (E N ,p N , 0,0), P H = P N -P L ={E H ,p H ) 

P L = ((E° L E N +Pop N cos 8 )/m Nl , (p E N cos9 + E° L p N ) /m Nl ,p sin 9 ,0) = (E L ,p L ) 
K = (uj, k cos 9, k sin 9 cos (p, k sin 9 sin ip) = (uj, k) 



where E° L 



[m 



m, 



m H ) /2m Nl and p 



m L are the L energy and 



momentum with respect to the N rest frame and 9q is a decay angle in the same frame. 
Now we compute the trace part. Using four-component Majorana spinors 14 , we get 



T 



polarizations 



X 



II L 



1-75 



Cv 



1 - 75 



i[#-f H + m N \C)x 



1-75 



1+75 



where C is the charge conjugation matrix with the properties 

Cv = u, {C, 7 n = 0, [C, 7 5 ]=0. 

Then eq. ()116j) becomes 

1+75 



;ii6) 



;ii7) 



T{K) 



-i Tr 



m Nl - 



'-P L m Nj [(l + a)#+blp] 



-2im Nl m N . [(l + a)P L -K + b E L ] 



:n8) 



Note that the term bJ/J is not put to zero by chirality projectors unlike a usual mass term. 

We have to perform the integral in d A K in eq. (|113jl . The most convenient technique is 
to use polar coordinates 9 and ip and integrate first in d cos 9 using the 5[(Pn — K) 2 — m 2 H ] 
and then integrate in dk using the 6[dL{v,k)]. As we discussed in section ET3| we can 
approximate the L dispersion relation with uj 2 — k 2 = m 2 L , finding the following solution 



cos 9 



m 



H 



m, 



m 



+ 2E N cu 



2kp 



k 



N 



lo £ — m 



L ■ 



(119) 



Imposing that | cos^l < 1 we find that uj must be comprised between two positive values 
co> m in and cj max . This implies that in the arguments of the ^-functions in eq. ()113j) we have 
u > and E N -uo > 0, so that expanding D%{P N - K)D%(K) + Djj(P N - K)Dl{K) gives 



N = [l + f B (E N -uj)][l-f F (uj)]-f B (E N -uj)f F (uj) 
= 1 + f B (E N -u)- f F {uj) - 2f B {E N - uj)f F {uj) . 



(120) 



14 The direction of the arrows for Majorana spinors in fig. 1161 is arbitrary, and one is free to choose it as 
a matter of convenience; the particular choice made dictates which Fcynman rules are used. 
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Figure 17: Iw\IqI-\\ in H — ► LN decay. 



At this point we are left with integrals in du and dtp. Adding also the contribution of the 
"wave" diagram the result is given by 

Im[(yt*%.] [ 2 * d^) /""»» J km Nl m N: 



(i2r 



27r(yty)u W y ^ \ds L /dk\ PN 

where 



|/ | 2 = 2P N ■ P L = m 2 Ni -m 2 H + m 2 L . (122) 
Py and Pw arise from iV} propagators in vertex and wave diagrams respectively 



1 _ 1 



mjj + ml- mjf - 2P H ■ K m z Ni - m% 



(123) 



appears in eq. (j!21j) multiplied by a factor 2 since in the wave diagram also the charged 
states of Higgs and lepton fields can propagate [28j . Finally, the relevant quantity for the 
Boltzmann equations is the average of €a over the phase space and the thermal distribution 
of Nr. 

P° dE N p N f B {E N ) f dcos9 e A (E N ,9 ) 
^ (T) = ,00 ,„ f ( „ , • (124) 

2 Jm. Nl dE NPN f B {E N ) 

The result is presented in fig. |B] In the limit m Nj ^> m Nl one has Py ~ Pw ~ — 1/™%. 
and the integrals in p and 9$ can be done analytically. Approximating /b(-EW) ~ e~ EN ^ T 
the explicit result is 



«*(0) - L ^ L ^ 1(1 + ^ ~ ^ + ^ + 2bEN] TK l{mN JT)\ dSL /8k\ ■ 



Computation of for N± at rest 

We have then performed the computation for the Ni decay in the simplified situation with 
N\ at rest in the plasma. In this case, the 5((Pn — K) 2 — m 2 H ) factor does not contain the 
angle 9 in the scalar product Pn ■ K. So we can proceed by integrating in dk and in duj 
using the two 5-functions. Moreover there is no dependence on the angles <p and 9q. The 
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integrals in dip and dcosOo are trivial and the integral in dcosO can be done analytically. 
In this way one obtains the result of section 13.31 The difference with respect to our full 
computation is not fully negligible, as shown in fig. 



Computation of e B 

Finally, the same technique of computation can be applied to the case of H — > LN%. There 
are, though, remarkable differences between the two decays. In H decay (see fig. IT7|) the 
particles in the loop could never go on-shell with the usual T = Feynman rules. This is 
due to the presence of the ^-functions in the cut propagators. However, using the finite-T 
rules we know that particles with negative energy can be absorbed from the bath so the 
process can have a CP asymmetry. We can make the cuts in fig. El where again cuts 
of very heavy particles have been neglected. The 5-functions have solutions identical to 
eq. ([119|) . in terms of the external momenta (see the notation of fig. [TTjl . However the 
masses are such that oo m i n < uJ mSuX < 0. So imposing | cos#| < 1, u has to lie between two 
negative values (it is exactly for this reason that with T = Feynman rules e B would be 
zero), while E N — u is positive. Then the products of ^-functions are different from the 
previous case 



Indeed eq. ()125j) is equal to (and not to 1) at T = 0. 

The result in this case has the same form of eq. (jl21)) . As already explained, since the 
final effect of e B in our scenario is very small, we computed it only taking the incoming H 
at rest, in order to simplify the computations. However, in this case it is not possible to 
give an analytic result as opposed to the case of N\ decay at rest. The reason is that this 
time the S((Pn — K) 2 — m 2 H ) factor contains always the dependence on the angle 8, since 
Ni is the outgoing particle and P/v is not of the form P N = (m Nl , 0, 0, 0). 

The most convenient order of integration is the same we followed in the Ni decay. We 
obtain again a result in an implicit form (even if there is no dependence on a and ip) 



D + H (P N - K)D+{K) + D H (P N - K)D- L {K) ex 

oc [f B (\E N - u \)][l - MM)] " [1 + fs(\E N - w|)]/f(M) 

= f B (\E N - u\) - MM) - 2f B (\E N - ^|)MM)- 



(125) 




x[f B (\E N - cu\) - MM) - Vb{\E n - u\)f F (\u\)] [P v (u) + 2P w (u)} , (126) 



where 5l is given by eq. (J115J) and I Q is the tree level rate for this decay 





m 2 Ni -m\. 



(127) 



The on-shell conditions fix 





H 



2E n lu 
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Finally Pv{&) and Pw(u) in this case are given by 

p y^) = — r^r — ^ — r- 5 v ' (128) 

Pn {m H + tul - m Nj z - 2m H u) 

Pw(cu) = * . (129) 

Pn {rriNi ~ m N 3 l ) 

Note that the tree level rate for this process in eq. ()127j) is small at high T as long as the 
value of rriH is near to mi,. For this reason e# is bigger than e^. In fact if |/o| 2 were 
C(1)T 2 at high T, then e# would be as big as ejvi- Instead, it goes like cT 2 , where c is 
a small number: c = (m 2 H — m 2 L )/T 2 , and so e# becomes c _1 times bigger than e^. In 
particular for the SM thermal masses c _1 is about 15. 



D Thermal correction to MSSM CP-asymmetries 

In the MSSM the situation becomes more complicated than in the SM because we must 
study both JVi and N\ decays, each having two possible decay channels, with each channel 
having more diagrams. Including thermal effects, the relevant masses satisfy 

m Nl23 —mfr , % = — 2~ - "7p a L = — — - 2a L . 

m N 1 z m N 1 

As previously discussed we use approximate dispersion relations for the fermions, and 
assume m N . ^> m Nl . For simplicity, in the MSSM case we compute the decay rates and 
the CP-asymmetries neglecting the velocity of the decaying particle with respect to the 
plasma. 

The CP-asymmetries in N\ and N\ decays are 



-N- 



V{Nx -> EL) - r(A>x -> EL) + r(A>! -> EL) - Y{N X -> EL) 

r(A>i -> hl) + r(A>i -> el) + r(iVi -> #l) + r(JVi -> #Z) 

i2 r (iVi -> EL)R t (N 1 -> iJL) + i? r (iVi -> EL)^^ -> #L) 



^(T = 0): 



i2 r (iVi -> #L) + i2 r (ATi -> i2X) 



e^(T) = e^(T = 0)R e (Ni — > = ejVj(T = 0)R e (Ni — > ifL) 
where 



r(i -> / at T = 0) v J y e(i -> / at T = 0) ' 

We have used the fact that at T = the two decay channels have equal zero-temperature 
width and CP-asymmetries ejvi(T = 0) = e N (T = 0) given in eq. (jUJ). 

We now give the explicit expressions for the R e . Consider first the Ni decays. Two 
more one loop diagrams (not plotted) contribute to the CP-asymmetry: a "vertex" and 
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a "wave" diagram with sparticles in the loop. The SM function R e (Ni — > HL) has been 
given in section l3~3l Using analogous self-explanatory notations we find, in the MSSM 

R^N^HL) = %J*L[ UL (l + a L ) + b L ][l + f H - f L -2f H f L \ 

m Nl 

k 2 

8— — ke-U + %) + b 6 ][l +n L -n 6 - 1n~ L n k \ 
m Nl 

where 

5 H = (m Nl - uo L f -k\-m\ 5 L = [(1 + a L )u L + b L } 2 - (1 + a L ) 2 k\ 
5 1 = (m Nl - Uh) 2 -k\-m\ 5 A = [(1 + a 6 )u 6 + b^] 2 - (1 + a k ) 2 k\. 

Re{N\ — > ifL) has to be evaluated at the values of ki and ujl which solve 5 B = 8l = 0, 
approximatively given by eq. (JHU), and at the values of /c^ and uig which solve 5jj = &i = 0, 
approximatively given by cj^ = (m^ + m 2 ^ — m|)/2m A r 1 and fc^ = (u^ — m 2 ^) 1 / 2 . 

Finally R e (Ni — > HL) = R e (Ni — > HL) since the CP-asymmetry in the two N\ decay 
modes is due to loops with the same virtual particles. 

For Ni decays, the situation is different. When Ni decays into two fermions (Ni — > HL) 
the imaginary part is obtained cutting two internal bosons, H and L. The decay rate 
is suppressed by two Pauli-blo eking factors, but its CP-asymmetry is enhanced by two 
stimulated-emission factors. 

R e (N, - HL) = 2-^-[l + f B (cu L ) + f B (u H ) + 2f B (u L )f B (u H )] (130) 

m N, 

where loh = { m % + m2 H ~ m |) l 2 ' m N 1 = m N 1 ~ w l an d k B = (oj 2 h — m 2 ^) 1 / 2 . 

When N\ decays into two bosons (N\ — > HL) the imaginary part is obtained cutting 
two internal fermions, H and L, and it is therefore given by a lengthy expression. The 
CP-asymmetry is suppressed by two Pauli-blocking factors, but the decay rate is enhanced 
by two stimulated-emission factors. 

k 2 r i 

R e (Ni — > HL) = 16— + a L )u + b L ][(l + a^Ug + b 6 ] + k 2 {l + a L )(l + a^)jx 



x [1 - f F {m^ - u) - /p(a;) + 2f F (m^ - u)f F (u)] 

where (u, k) is the quadri-momentum of L, ug — irift — oj is the energy of if and 
^ = [(1 + a^)(m^ -u) + bg} 2 - (1 + a^) 2 /c 2 <5 L = [(1 + a L )uj + b L ] 2 - (1 + a L ) 2 A; 2 . 

The expression should be evaluated at the values of oj and k which solve 5l = Sg = 0. 
Using the approximate on-shell condition, they are given by 

u = ( m yv! + m l ~ m \)/ 2m Nv k={u 2 - m 2 L ) 1/2 . 
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d5 H /doj L dS L /doj L 
dS H /dk L d5 L /dk L 



dS L /du s dS 6 /du 6 
dS- L /dkH dS 6 /dk 6 



d5 L /du dSfj/doj 
d5 L /dk dS^/dk 



'131) 



Finally, the thermal corrections to the decay rates are given by 



RriN^HL) = (l + f B (E H ))(l-f F (E L ))\^ 2 (l,a H ,a L )(l-a H + a L ) 

Rt(Ni — > HL) = (l-/ F (^))(l + / B ( J E; z ))A 1 / 2 (l,a^,a L )(l + a i? + a z ) 

Rv{Ni — > HL) = (l-/ F (^))(l-/ F (E L ))A 1 / 2 (l,a^,a L )(l-a^-a L ) 

Rt{N\ — > HL) = (l + / B (^))(l + / B (^z))A 1/2 (l,^,a z ). 

The 1 ± Jb,f factors take into account Pauli blocking or stimulated emission (the thermal 
distributions must be evaluated at the energies of final state particles) while the other 
factors arise from thermal corrections to kinematics. 



Note added Some of our preliminary results appeared in [52] • However, in that paper 
sneutrino reheating was not correctly included (we also take into account thermal correc- 
tions and proper subtraction of on-shell scatterings). 
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